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Abstract 



The prevailing theory for the origin of cosmic magnetic fields is that they have been 
amplified from insignificant seed values to their present values by the turbulent dy- 
namo inductive action driven by the plasma turbulent motions in the protogalactic 
and galactic medium. Up to now, in calculation of the turbulent dynamo, it has been 
customary to assume that there is no the back reaction of the magnetic field on the 
turbulence, as long as the magnetic energy is less than the turbulent kinetic energy. 
This assumption leads to the kinematic dynamo theory that has been well developed 
in the past. 

However, the applicability of the kinematic dynamo theory to protogalaxies is 
rather limited. The reason is that in protogalaxies the temperature is very high, 
and the viscosity is dominated by ions. As the magnetic field strength grows in 
time because of the dynamo action, the ion cyclotron time becomes shorter than 
the ion collision time, and the plasma becomes strongly magnetized. As a result, 
the ion viscosity becomes the Braginskii viscosity, and the magnetic field starts to 
strongly affect the turbulent motions on the viscous scales. Thus, in protogalaxies 
the back reaction sets in much earlier, at field strengths much lower than those which 
correspond to energy equipartition between the field and the turbulence, and the 
turbulent dynamo becomes what we call the magnetized turbulent dynamo. 

The main purpose of this thesis is to lay the theoretical groundwork for the mag- 
netized turbulent dynamo. In particular, we predict that the magnetic energy growth 
rate in the magnetized dynamo theory is up to ten time larger than that in the kine- 
matic dynamo theory, and this could lead to the dynamo creation of cluster fields. We 
also briefly discuss how the Braginskii viscosity can aid the development of the inverse 
cascade of magnetic energy, which happens after the energy equipartition time. 
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Chapter 1 



Introduction 



1.1 Galactic and Extragalactic Magnetic Fields: 
Observational Results 

One of the most important and challenging questions in astrophysics is the origin of 
galactic and exragalactic magnetic fields. The existence of interstellar magnetic fields 
in our Galaxy was first proposed by Alfven in 1943 In 1949 the polarization of 
starlight was observed by Hiltner and independently by Hall and Mikesell [|17| . 



The starlight polarization was interpreted as a consequence of light scattering by in- 
terstellar dust grains aligned in the galactic magnetic fields. In the same year Fermi 
employed galactic magnetic field for the acceleration of the cosmic rays and for their 
confinement in the Galaxy [|1^. In 60's the Faraday rotation and the Zeeman ef- 



fect were measured for different sources distributed over the sky |T8|, In 1970 
Mathewson and Ford found a large-scale magnetic field in Magellanic Clouds by mea- 
suring the polarization of stars in the Magellanic System In 1974 Manchester 
first correctly measured the galactic magnetic field in the vicinity of the Sun (within 
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~ 2 Kpc) by measuring the Faraday rotation for thirty eight nearby pulsars . He 
found that the local field has a longitudinal component in the Galactic plane. In 
1989 Rand and Kulkarni confirmed and significantly improved the measurements of 
the local magnetic field They analyzed the Faraday rotation measures for nearly 
two hundred pulsars within ^ 3 Kpc, and they found that the uniform component of 
the local magnetic field has a strength of ~ 1.6 fiG towards a galactic longitude of 96°, 
with a reversal of the field at a distance about 0.6 Kpc towards the galactic center. 
The random field component was estimated by Rand and Kulkarni as ~ 5 fiG. Since 
early 90's a significant progress has been made in measuring magnetic fields in our 
Galaxy and in other spiral galaxies, by measuring polarization of starlight, by mea- 
suring the Faraday rotation from pulsars and extragalactic radio sources, by detecting 
synchrotron radiation from relativistic electrons, and by measuring the Zeeman split- 
ting of spectral lines; see recent excellent review papers 0, ^ The observations 
indicate that galaxies possess magnetic fields with strengths of several microgauss (up 
to several tens of microgauss). These fields have uniform components, with strengths 
comparable to those of the random components (e. g. in the vicinity of the Sun the 
uniform/random field strength ratio is about 1/2). The uniform field components 
lie in the galactic disks and have typical correlation lengths from several hundreds 
parsecs to one kiloparsec. Figure ^]l] shows the magnetic field pattern for the galaxy 
M51. However, most galaxies do not show such a nice regular field pattern as M51 
does. As a result, the question about the prevailing geometrical structure of magnetic 
fields in spiral galaxies is still open, see Figure |1.2| . The magnetic field in Milky Way 
exhibits reversals with radius, which suggests a bisymmetric spiral structure for it, 
see Figure |T]2|(B). 

It is very important to know that magnetic fields may exist in very young galaxies 
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M51 6cm Total Intensity + Magnetic Field 




Figure 1.1: Magnetic field vectors (white bars) for tlie galaxy M51 (measured by 
fitting a model of the polarized synchrotron emission to radio observations of M51). 
The map is superimposed on a M51 image in 6 cm. The magnetic field lines approxi- 
mately follow the spiral structure. By courtesy of Beck, Horellou, Neininger, and the 
MPIfR collaboration [H. 
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A B 

azimuthal azimuthal 




Figure 1.2: Two possible magnetic field configurations in disk galaxies, face-on view. 
The bold lines are magnetic field lines with the field direction shown by the arrows. 
A: axisymmetric spiral structure, no radial reversals of the field (this structure is 
predicted by the galactic dynamo theory); B: bisymmetric spiral structure, there are 
radial reversals of the field (this structure arises in the rotating galactic disk if the 
field has primordial origin |]58|). 
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(at high redshifts) [Q. In 1992 Wolfe, Lanzetta and Oren found that the probabihty 
for Faraday rotation for radio QSOs (quasi-stellar objects) is significantly higher in 
damped Lya systems [^. These systems are associated with early forming galactic 
disks. Wolfe et al. estimated the magnetic field in two damped Lya systems with 
few microgauss. More recent extensive observations have confirmed this 
result. Thus, magnetic fields with similar spatial scale and strength to those in the 



local universe exist at redshifts above 2 = 1, see review papers ||2^, |3^, The 
field strength estimates all fall in the range of 1 to 5 /iG, and there is an absence 
of observational evidence for cosmological evolution of magnetic field strengths! In 
addition, there also exist some recent observational data which indirectly indicate the 
existence of rather strong magnetic fields in the distant past. The observations of ^Li 
and ''Li abundances in old metal-pure halo stars by Lemoine et al. indicate a massive 
production of lithium isotopes by cosmic rays in a very early phase of the Milky Way 
Galaxy ||30| . Without magnetic fields the cosmic rays would escape and would not be 
able to produce the lithium isotopes. Other indirect evidence for magnetic fields in the 
past is the primordial star formation process, which must have happened according 
to the observed metallicities of globular clusters. To form stars at early times (z ~ 3), 
primordial magnetic fields are believed to be required in order to remove significant 



angular momentum of self-gravitating gas by the magnetic braking effect [40 



Strong magnetic fields have been observed not only in galaxies, local and at high 



redshifts, but also in clusters of galaxies |I4|, El, |58[ . The most recent direct observa- 



tions of the Coma Cluster by Fusco-Femiano et al. [|T6[ manifest a volume averaged 
intracluster magnetic field of ~ 0.15 /xG in the cluster. Fusco-Femiano et al. used the 
inverse Compton scattering of relativistic electrons on the cosmic microwave back- 
ground photons to model the hard nonthermal X-ray radiation flux from the Coma 
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Cluster. As a result, they estimated the electron density in the cluster. Then, they 
used the synchrotron radio flux to estimate the fleld strength. In 2000 Sarazin and 
Kempner revised the results of Fusco-Femiano et al. . Sarazin and Kempner used non- 
thermal bremsstrahlung radiation models for the hard X-ray emission from the Coma 
Cluster, and depending on the model they used, they estimated the fleld strength 
to be ranging from ~ 0.4 yuG (as also implied by equipartition in the radio halo) to 
~ 6 jiG (as the typical fleld in individual galaxies in the cluster) [^] . 

Thus, galaxies, early galaxies (at z ~ 2 redshift) and galaxy clusters do have strong 
magnetic flelds, with the strengths 0.1 /iC < i? < 10 /iC! Where did these strong flelds 
come from? 



1.2 The Origin of Galactic and Extragalactic 
Magnetic Fields: Primordial and 
Galactic Dynamos 

The prevailing theory for the origin of strong cosmic magnetic flelds is that they were 
produced by the turbulent dynamo inductive action driven by the fluid motions in 
galactic and/or protogalactic medium. The turbulent dynamo works as follows. The 



astrophysical plasmas are very hot and they have low densities (e. g. see Table |LT 
Therefore, the effect of electrical resistivity is negligible over an extremely broad range 
of scales. As a result, the magnetic fleld lines are frozen into plasma [^, Because 
in a turbulent plasma the distance between any two inflnitesimally neighboring points 
increases exponentially in time (the Kolmogorov-Lyapunov exponentiation), the mag- 



netic fleld lines are stretched and the fleld strength grows exponentially fast [22, 55 



see the left plot in Figure 1.3. There is an additional mechanism for the growth of 



the mean field strength in the case of three-dimensional turbulent dynamos, which 



is called the Zeldovich "figure 8" mechanism p9l pH]. This mechanism is shown in 



the right plot of Figure |1.3| , where the closed magnetic field lines are frozen into the 
torus. The field strength is increased exponentially fast in time by repetition of first, 
stretching the torus, second, twisting it into "figure 8", and third, folding it. 

Yet full understanding of all stages of production of the strong galactic and ex- 
tragalactic magnetic fields by the turbulent dynamo action has not been achieved. 
There are two alternative theories on how and when the fields have been produced. 

The first theory, the galactic dynamo theory, also known as the a-Q dynamo the- 
ory, states that the fields have been primarily amplified in differentially rotating galac- 
tic disks after the galaxies had been formed. This theory considers the evolution of a 
large-scale mean magnetic field in a galactic disk. The turbulent dynamo action and 
the differential rotation of the disk are parametrized by transport coefficients. These 
transport coefficients enter the evolution equation for the mean field and they repre- 



sent the destruction and the induction of the mean field 0, ^ |3g, ^7} |5y, ^ |5J, p5 
The typical time scale of the mean field growth in the galactic disk is approxi- 
mately equal to the disk rotation period, ~ 300 million years. The galactic dy- 
namo theory is a beautiful theory, but it involves several crucial unsolved prob- 
lems [|^, |19|, |5^. First, in the a-Q theory it seems to be extremely 
difficult to expel a fraction of the magnetic flux from the galactic disk in order to 
produce the net magnetic flux [^, The second problem is that the effect of 
small-scale fields on the amplification of the mean magnetic field remains unclear. 
These small-scale fields are amplified on the smallest turbulent eddy turnover time 
scale, much faster than the mean field is being amplified (see Figure |1.5| ). Thus, the 
small-scale fields must saturate and possibly significantly alter the a-Q dynamo ac- 
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c* - 'd 

y stretch 




Figure 1.3: Left plot: stretching of magnetic field lines in a turbulent plasma because 
of the Kolmogorov-Lyapunov exponential divergence of closely neighboring points. 
Right plot: three-dimensional Zeldovich "figure 8" dynamo mechanism. 
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tion well before the mean field is amplified appreciably The third problem is 



that the galactic dynamo theory requires rather large seed fields in the distant past, in 
order to successfully build up the fields observed at the present time. The Biermann 
battery itself and other similar effects do not seem to be able to provide the required 
seed fields Finally, as we wrote in Section |1 . 1| , the observational data indicate 



that galaxy clusters and galaxies at low and high redshifts possess magnetic fields 
of up to several microgauss. It is hard to explain these strong fields by the galactic 
dynamo theory [Q. On one hand, the galactic dynamo just did not have enough 



time to build microgauss fields in early galaxies at, say, redshift z ~ 2. On the other 
hand, the fields in clusters of galaxies are several orders of magnitude larger than 
the fields that would be obtained from galactic ejecta alone. Thus, sufficiently strong 
initial magnetic fields must be generated during the pregalactic era! 

In this thesis we accept the second theory for the origin of cosmic magnetic fields, 
the primordial dynamo theory, which states that the galactic and extragalactic mag- 
netic fields have primarily been produced in protogalaxies, i. e. before the galaxies 
were formed [^, ^ ^ Of course, these fields were subsequently modified in the 
rotating galactic disks after the galaxies were formed f^. 

In order to understand how the magnetic fields can be built up in protogalaxies, 
let us first discuss the physical conditions that were present there. The typical values 



of physical parameters in a protogalaxy are given in Table |L^, and the hierarchy of 
scales is displayed on the x-axis of Figure |1.4| . The most important facts are the 
following. 

• First, the gas is very hot in a protogalaxy, so it is fully ionized. Therefore, 
the viscosity is dominated by ions, not by neutrals. In addition, the Spitzer 



resistivity |^ is tiny, and thus, the resistive cutoff scale for the magnetic field. 
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Parameter 


Notation 


Value* 


Scahng* 


Physical Quantities 


total mass, g 


IVl 


Ji X iU 


^ , 1 nl2 TV /r 
~ lU J-VIq 


total/baryon mass ratio 


t 
? 


iU 




temperature, K 


T ~ ~ Te 


2 X 10^ 


ML-^ 


ion & e~~ density, cm~^ 


n 


5 X 10-* 




neutral density, cm~^ 









ion thermal speed, cm/s 


Vt 


2 X 10^ 




ion viscosity, cm^/s 


= z/j ~ V^ti 


5 X 10^6 




thermal speed, cm/s 


{mi/meY''^VT 


9 X 10^ 




viscosity, cm^/s 




1025 


^^3/2^1/2 


neutral viscosity, cm^/s 









magnetic diffusivity, cm^/s 


Vs 


8 X 10* 




smallest eddy speed, cm/s 




2 X 10^ 


^1/4^3/4^-1/2 


Dimensionless numbers 


hydrodynamic Reynolds 


R ~ VtL I V 


3 X 10* 




magnetic Reynolds 


Rm ~ VTL/r]s 


2 X 10^^ 




Prandtl number 


Pr ~ v /lis 


6 X 10^1 




field/smallest eddy energy 


B-/ inni^^nVj 


3 X l{f^B''- 


e'-M-''-L'B^ 


Length Scales 


system size, cm 


L 


6 X 10^3 


~ 0.2 Mpc 


viscous cutoff scale, cm 




3 X 1020 


e3/*M3/*L 


ion mean free path, cm 


= VtU ~ R-^L 


7 X 10^9 




ion gyroradius, cm 


Ti = Vr/LJi 


2 X W/B 


jyji/2^-1/25-1 


e~ mean free path, cm 


~ {me/rriiY/'^Xi 


2 X 10^^ 




e~ gyroradius, cm 




50/5 


^1/2^-1/2^-1 


resistive cutoff scale, cm 


27rfc->Pr-V2i?-3/42^ 


4 X 10^ 


^1/4^-3/4^^3/2 


Debye length, cm 


(fcBT/27me2)i/2 


6 X 10^ 




Time Scales 


gravitational collapse time, s 


~ L/Vx 


3 X lO^*^ 




largest eddy turnover time, s 


~ L/Vt 


3 X 10^^ 




smallest eddy turnover time, s 


~ R-y^L/Vr 


2 X 10^* 


^1/2 ^3/2 


ion collision time, s 


ti 


3 X 10^2 


^^1/2^3/2 


ion cyclotron frequency, s~^ 


uji = eB/rriiC 


9 X WB 


B 


e~ collision time, s 




7 X 10^° 




e~ cyclotron frequency, s~^ 


{mi/me)oJi 


2 X WB 


B 



Table 1.1: Physical parameters in a protogalaxy calculated for a fully ionized hy- 
drogen plasma, assuming definite values for the parameters printed in boldface. * The 
field B in column "Value" is expressed in gauss, the Coulomb logarithm is taken to 
be 30. 
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Figure 1.4: The Kolmogorov spectrum of kinetic energy, 1^, is shown by the thick 
sohd hne. The x-axis has the wave number k plotted on it, and demonstrates the 
hierarchy of scales in a protogalaxy: L is the system size, 27rA;~^ is the viscous cut- 
off scale, Aj is the ion mean free path, and 27rfc~^ is the resistive scale. Note that 
L ^ 27r/c~^ ^ Aj ^ 27ik~\ The equations for a nonresistive incompressible MHD 
turbulence are valid only inside the range of scales shown by the shaded area (see 
also footnote Q on page |12|) . 
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271 k~^, is extremely small compared to the viscous cutoff scale for the turbulent 
velocities, 27rA;j7^. As a result, we can neglect resistivity as long as we consider 
scales ^ 2Tik~^, see Figure |L4 . 



Second, the Reynolds number is very large. Therefore, the turbulence develops 
over a very broad range of scales, starting at the largest scale, L, and ending 
at the viscous cutoff scale, 2T[k^^ (see Figure |1.4|) . The velocities at the largest 



scale are of the order of the sound speed, which is approximately equal to the 
thermal speed. All turbulent velocities at smaller scales are smaller. As a result, 
we can treat the plasma as incompressible on scales 2T{k~^ <^ L. 

Third, the ion mean free path Aj is considerably shorter than the viscous cutoff 
scale for the turbulent velocities, 27: k~^. The Debye length and the resistivity 
scale are even much shorter than Aj. Thus, we can use the single- fluid magneto- 
hydrodynamic (MHD) equations for description of a nonresistive incompressible 
plasma on scales Aj < 27r/c~^ ^ Q 

Finally, it is very important that the cyclotron period of ions in the magnetic 
field is shorter than the ion collision time, lu^^'^ <^ ti, provided that the field 
strength is larger than ~ 10^^^ G (see Table |0|). On the other hand, the 
energy of the magnetic field becomes comparable to the kinetic energy of the 
smallest turbulent eddies (which are on on the viscous cutoff scale) if the field 
exceeds ~ 10~^G. As a result, there is a very broad range of magnetic field 
strength values at which the magnetic pressure and tension are still negligible, 
while the presence of the field is already important. This is because the plasma 



^ We should consider scales larger than Ai in order to use MHD description for the plasma. 
However, for estimation purposes, one can use MHD equations even for smaller scales, 27rfc^^ ^ A^, 
if one reduces the molecular viscosity by the ratio of the scale to the ion mean free path, i. e. by 
factor l/{kXi) < 1. 
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is strongly magnetized, and the magnetic field controls the microscopic motions 
of ions, so that the plasma viscosity is different from that in a field-free plasma. 

The primordial dynamo theory is as follows. We believe that there are five major 
stages of the production of the strong magnetic fields. Four of them happen in 



protogalaxies, and one happens in galaxies, see Figure |1.5| . These stages are the 
following: 

1. During the first stage, in a protogalaxy undergoing gravitational collapse, the 
Biermann battery builds a seed magnetic field linearly in time. The Biermann 
battery requires the pressure to be non-barotropic, which is achieved behind 
hydrodynamic shocks driven by the gravitational instability The resulting 
seed field is of the order of IO^^Liq-is G at the largest scales [0, ^, |D|. 



The Biermann battery inductive action is proportional to the vorticity of the 
turbulence, which is R^/^ times larger at the viscous scales compared to the 
vorticity at the largest scales (i? 1 is the Reynolds number). Therefore, the 
seed field at the smallest scales of the turbulence (the viscosity scales), is about 
~ 10^^^ G. The time scale of the Biermann battery action is approximately 
equal to the free-fall time in a protogalaxy, ~ 1 billion years. 

2. During the second stage, when the plasma is unmagnetized and uj^^ ^ tj, the 
seed field is exponentially amplified by the kinematic turbulent dynamo inductive 
action. The kinematic dynamo builds the field up to approximately 10^^^ G, 
when the plasma becomes magnetized, uj^^ <^ ti, and the field becomes strong 
enough to affect microscopic processes in plasma, such as the viscosity. The time 
scale of the kinematic dynamo is very short, it is of the order of the smallest 
eddy turnover time, ~ 5 million years. The main assumption of the kinematic 
dynamo theory and, basically, the strict definition of it, are that the growing 
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Figure 1.5: There are five major stages of the production of strong galactic and ex- 
tragalactic magnetic fields, as the field strength grows from zero up to < 10^^ gauss. 
The results of this thesis apply to the magnetized turbulent dynamo stage in proto- 
galaxies, shaded on the plot. This stage makes the largest contribution to the built 
up of the magnetic fields. 
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magnetic field stays so weak, that it does not affect the fluid motions, i. e. that 
there is no the hack reaction of the field on the turbulence. The kinematic 
dynamo has been intensively studied during the last several decades (see the 
references given in Section p.lD . Up to now it has been customary to assume that 
there is no the back reaction, and the main assumption of the kinematic dynamo 
is valid, as long as the magnetic energy is less than the turbulent kinetic energy, 
so that the Lorentz forces are small. However, the plasma becomes strongly 
magnetized and the back reaction sets in much earlier, at field strengths much 
lower than those which correspond to the energy equipartition between the field 
and the turbulence (see Figure pTS] ). Thus, the applicability of the kinematic 
dynamo theory to protogalaxies is rather limited. In Section we discuss the 
main results obtained for the kinematic dynamo theory. The results we obtain 
in this thesis reduce to those of the kinematic dynamo theory in the limit of 
very weak field strength (i. e. when the plasma is unmagnetized, and there is 
no the back reaction). 

3. The third stage starts when the field grows above ~ 10^^^ G, the ion cyclotron 
time in the magnetic field becomes shorter than the ion collision time, and the 
plasma becomes strongly magnetized. As a result, the magnetic field starts to 
strongly affect the dynamics of the turbulent motions on the viscous scales 0, by 
completely changing the viscosity (see Section p^^), despite the fact that the field 
energy is still very small compared to the kinetic energy of the turbulence! We 
call this stage as the magnetized turbulent dynamo (see Figure |L^). In previous 

theories this stage has not been recognized. The main goal of this thesis is to 

^ Note that the turbulent motions on the viscous scales are the most important in the dynamo 
theory. This is because the magnetic field is primarily amplified by the turbulent eddies on these 
scales. These eddies are the smallest ones, they have the shortest turnover times and produce the 
largest velocity shearings ]27[ |. 
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construct a theoretical model for it. The time scale of the magnetized dynamo is 
roughly about the same as the time scale of the kinematic dynamo, ~ 5 million 
years. 

4. So far the field scale is of the order of the viscous scale or less, and the magnetic 
field is incoherent in space. The fourth stage starts when the magnetic field 
strength grows up to ~ 10^^ G. The field energy becomes comparable to the 
kinetic energy of the smallest turbulent eddies (the energy equipartition), and 
the Lorentz forces become dynamically important in the plasma. During this 
stage the turbulent motions are dissipated by the growing field, the turbulent 
energy spectrum becomes truncated at larger and larger scales, and the turbu- 
lent energy is eventually transfered into a large-scale strong magnetic field with 
its energy comparable to the kinetic energy of the fluid motions on the largest 
scales in the protogalaxy (which is approximately the same as the thermal en- 
ergy). This process is called the inverse cascade g, |, |7l |3|, |ll • We discuss it 



only qualitatively in Chapter ^ The time scale available for the inverse cascade 
process is of the order of the largest eddy turnover time, ~ 1 billion years. The 
inverse cascade may not have time to amplify the field up to microgauss values, 
which are observed in galaxies. A crucial question is how far it goes. This thesis 
addresses this question and is concerned with the rate of the magnetic field built 
up by the magnetized dynamo. 

5. Finally, the fifth stage is the galactic dynamo, which happens in the differentially 
rotating galactic disc after the galaxy is formed (see the discussion of the galactic 
dynamo above and Figure |1.5| ). This process modifies the strong field that was 



built up in the protogalaxy in advance The galactic dynamo may also 



further amplify the mean magnetic field in the disk by a factor of three or so, 
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on a time scale of the order of the rotation time in the galaxy, ~ 300 million 
years. The galactic dynamo theory is beyond the scope of this thesis. 

As we discussed in the beginning of this section, the galactic mean field dynamo 
theory seems to be unable, by itself, to build magnetic fields from the seed fields of ~ 
10~^^ gauss, which are provided by the Biermann battery action, up to > 10~^ gauss, 
which exist in the present Universe. The observational evidences and the theoretical 
problems of the galactic dynamo theory, discussed above, persuade us to accept the 
primordial dynamo theory, which states that the observed cosmic fields were primarily 
built up in protogalaxies. The primordial dynamo theory seems to be free of the 
shortcomings of the galactic dynamo theory. The turbulent dynamo in protogalaxies 
does seem to be able to amplify the seed field up to ~ 10^ G in a time interval 
< 10^ years, which is negligible in comparison to the Hubble time. However, up to 
now, there was one major deficiency in the primordial dynamo theory. It is well 
known that the turbulent dynamo builds the magnetic fields first on small subviscous 



scales [^, while the observed cosmic fields have rather large correlation lengths. 
Therefore, the magnetic field lines must be unwrapped on the small scales by the 
Lorentz tension forces, while the field energy is transfered and amplified on larger 
scales during the inverse cascade stage. Up to now, there was a strong argument 
against the possibility of such unwrapping. This argument is as the following. Both 



theoretical calculations and numerical simulations |Tl| , |^ show that the magnetic 



field, built up by the standard dynamo action, has folding structure with the wave 
number perpendicular to the field lines much larger than the wave number parallel to 
the field lines, as demonstrated in Figure |2TT|. Therefore, the Lorentz tension forces. 



which are proportional to the parallel wave number, are not very large. As a result, 
the standard isotropic viscous forces, calculated for an unmagnetized or partially 
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ionized plasma, can balance the Lorentz tension forces 0, and block the unwrapping 
of the field lines (see Chapter At the same time, the perpendicular wave number 
can grow in time, and the perpendicular scale of the field lines can unrestrictively 
decrease 0. 

This anti-unwrapping argument [jlT| applies only if the viscous forces are isotropic. 



However, in the case of the magnetized turbulent dynamo in a fully ionized plasma, 
the magnetic field controls the microscopic motions of particles on the viscous scales. 
The transport of ion momentum across the field lines is inhibited, and the viscosity 



stress is given by the Braginskii tensor [see equation ( |2.28| ) and the discussion before 
it]. As a result, in the magnetized dynamo theory the viscous forces do not prevent 
the unwrapping of the magnetic field lines by the Lorentz tension forces on the small 
scales and do allow the field to become converted into a large-scale field during the 
inverse cascade stage in a protogalaxy (see Chapter || for more discussion). On the 
other hand, the theory of turbulent dynamos with the Braginskii viscosity has not 
been worked out. It turns out that the dynamo theory is qualitatively modified 
by including the Braginskii viscosity. The main purpose of this thesis is to lay the 
theoretical groundwork for the magnetized turbulent dynamos. 

1.3 The Thesis Preview 

Let us briefly outline the thesis content. 

In Chapter H we give the main equations and the main results of the kinematic dy- 
namo theory (in Section |2.1| ), and we formulate the basic equations of the magnetized 

dynamo theory (in Section p.2|) . 

The magnetic pressure is balanced by the hydrodynamic pressure in an incompressible plasma. 
^ Until it reaches the scales at which resistivity becomes important, see Figure L4 on page |ll|. 
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In Chapter ^ we calculate the statistics of turbulent velocities, V, in strongly 



magnetized plasmas. In Section we write down the quasilinear expansion in time 
of the MHD equations for both the velocities and the magnetic field, similar to the 
expansion used by Kulsrud and Anderson P, ^ . In section |3.2| we make use of the 
Laplace transformation in time to find the turbulent velocities. We assume that the 
tensor 6q,/3, which is the product of two unit vectors in the direction of the magnetic 
field, see equation (|2.38|) , can be taken to be constant in space in the beginning of 
the expansion in time, ba/3 = habp = const at zero time. This is our first hypothesis. 
It basically relies on our assumption that in the case of the magnetized turbulent 
dynamo the magnetic field has a folding structure similar to the one that exists in 
the case of the kinematic turbulent dynamo (see Figure |2.1j and the discussion in 
Section p.2| ). We find that there are velocity modes which are not damped by the 
Braginskii viscous forces. These undamped velocity modes grow unrestrictively in 
time unless we incorporate the non-linear inertial terms of the MHD equations into 



our quasilinear expansion. In Section |3.3| we argue that the non-linear damping of 
the growing velocity modes may be included into our theory by allowing for rotation 
of velocity vectors relative to the magnetic field unit vectors |. This is our second 
hypothesis. In section [3l^ we again solve for the turbulent velocities by now including 
our effective rotational damping into the MHD equations and by making use of the 
Fourier transformation in time. As a result, we obtain the correlation functions for 
the turbulent velocities in a strongly magnetized plasma, similar to those given by 
equations ( |2.4| ) and ( p.5| ) for the Kolmogorov velocities in the kinematic dynamo 
case. We find that, contrary to the Kolmogorov velocities, the turbulent velocities 

in the magnetized plasma are strongly anisotropic, as one might expect, because the 

^ This rotation is essentially due to Coriolis forces that make the velocity rotate differently than 
the magnetic field direction. 
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magnetic field sets "a preferred axis in space" . 

In Cfiapter ^ we use tlie correlation functions for the turbulent velocities found 
in Chapter § to calculate the energy spectrum of random magnetic fields in the mag- 
netized plasma. We start with calculations of the total magnetic energy growth rate 
in Section |4.1| . In Section |4.2| we derive the mode coupling integro-differential equa- 
tion for the evolution of the magnetic energy spectrum. Our mode coupling equation 
is a more general version of equation ( p.21j ) of Kulsrud and Anderson obtained in 



the kinematic dynamo theory '^Tj. In Section |4.3| we consider the magnetic energy 
spectrum on small subviscous scales. On these scales our mode coupling equation 
greatly simplifies and becomes a homogeneous differential equation, similar to the 
equation ( p.25| ) of Kulsrud and Anderson. We find the Green's function solution of 
our differential equation. 

Finally, in Chapter]^ we give our conclusions. We discuss the possibilities of further 
research on the magnetized turbulent dynamos. We also discuss the peculiarities of 
the inverse cascade in a strongly magnetized turbulent plasma. 
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Chapter 2 



Basic Dynamo Equations 



2.1 The Kinematic Turbulent Dynamo 

The main assumption of the kinematic turbulent dynamo theory is that the magnetic 
field is weak and it does not affect the turbulent motions in plasma. Thus, the 
magnetic field B, frozen into plasma, is evolved as a passive vector by the turbulent 



velocities U 29 : 



dtBa = Ua,l3Bp — UpBa^p. (2.1) 



Here and below we always assume summation over repeated indices. In order to 
shorten notations, we frequently use dt == d/dt, and spatial derivatives are assumed 
to be taken with respect to all indices that are listed after " , " signs []. We also 
assume that the turbulence is incompressible (see Section |1.2| ). In the kinematic 
dynamo theory the turbulent velocities U are assumed to be hydrodynamic and to 

be independent of the magnetic field. 

^ For example, dtB^ = dBa/dt, {UaBp)^^ = B/j{dUa/dxj) + Ua{dBp/dx^), and Ua^p-y = 
d'^Ua/dxpdx-y. 
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A vast number of papers have been written about the kinematic turbulent dynamo 
in order to answer three main questions: first, "How fast is the total magnetic energy 
built up by the kinematic dynamo action", second, "What is the magnetic energy 
spectrum on different scales?" , and third, "What is the resulting geometrical structure 

, |ni n, n, 13, n, n, |9|, n, n, 0, ^ |55|. It 



of the magnetic field?' 



turns out that in most cases the evolution of magnetic field is entirely determined by 
the specified two-point statistics of the turbulent velocities. 

Following the assumptions usually made in the kinematic dynamo theory, we 
assume that the turbulence is incompressible, homogeneous, isotropic and stationary 
(as it actually is in protogalaxies) . We also neglect the helical part of the turbulence 0. 
We assume that the statistics of the Fourier coefficients of the turbulent velocities, 



L/2 



L3 



-L/2 



Ukait) e*-^* dt, 



(2.2) 
(2.3) 



is given by the following formulas [27|: 



(f/ka(^)) = 0, (2.4) 

{UU^)U^.'p{uJ')) = {U\a{-^)Ui.'p{oo')) = ^fc5i5k',-k5(^' + ^). (2.5) 



Here and below (...) means ensemble average over all realizations of the turbulence, 

^k',-k is the three-dimensional Kronecker symbol (equal to unity only if k' = — k), 

^ The helicity is negligible in galaxies on the scales of the smallest turbulent eddy, which is the 
principal driver of the field evolution |27j . The helicity in protogalaxies is even smaller than that in 
galaxies. 
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5{uj' + uo) is the Dirac (5-function, 



= ^o^p - kakp, (2.6) 



Sap is the one-dimensional Kronecker symbol, and k = k/k is the unit vector along 
the wave vector k. It is important that function J^^, which stands for the normal 
(non-helical) part of the turbulence, depends only on the absolute values of uj and k. 
Note that equations ( |2.2| ) and ( |2.3| ) represent the discrete three-dimensional Fourier 
transformation in space and the continuous one- dimensional Fourier transformation 
in time (see Appendix 0) , so that the inverse Fourier transformations are 

r 3 °° 

Ua{t,r) = E^ka(t)e*'^'-= (— ) J U^{t,k)e'^''d% (2.7) 

^ — oo 

1 °° 

U^4t) = ^ I U^aHe-'^^'du. (2.8) 



-oo 



Here the summation is done over discrete values of vector k: kx = {27r /L)nx, ky = 
(2,71 /L)ny, kz = {2TT/L)nz, & Z, Uy & Z, Uz G Z] and this summation can be re- 
placed by integration of the appropriately defined function Uait, k) that is continuous 
over k, see Appendix 0. Below we will use both summation and integration over k, 
whichever is more convenient to use. 

Let us further assume that the turbulence is Kolmogorov, and that the time 
correlation function of the turbulent velocities has an exponential profile i. e. that 

(U(t)U(t')) oc e-l*-*'l/", where r is the eddy decorrelation time which depends 
Using a Gaussian time cor relation profile, (U(i)U(t')) cx e-(*-*') , would be more appro- 



3 



priate. In this case equation ( 2.10 ) would become J^^k = Jak^'^ However, we prefer the 

exponential profile because it is easier to deal with. (In particular, for Gaussian integrals it is not 
possible to close integration contours at infinity in the complex plane.) 
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on k, 

Here, /cq = 27r/L is the smallest wave number of the turbulence, and Uq ^ Vt is the 
largest eddy velocity. In this case we have 



J.k = — (2.10) 

1 + T'^UJ^ 



iUo/2ko)ik/kor''^\ ko<k<k,, 
Jok ~ { (2-11) 
0, k < ko, k > ku. 



and, carrying out the inverse Fourier transformations of equation (p.5|) , in time and 
in space, we have 



(t/ka(t)t/k'MO) = ^e-l*-*'l/^5i/3 5k',-k, (2.12) 

zr 

(f/„(t,r)f/,(t',r')) = Ey^e-l*-'l/^5„^,e^'^(-'). (2.13) 

k 



The Kolmogorov kinetic energy spectrum, shown by the thick solid line in Figure 1.4, 



is 27 



, {Ullko){k/k,)-"\ ko<k<k,, 
I{k) ^ { (2.14) 

0, k < ko, k > k,y, 

so that the total kinetic energy of the fluid motions, per unit mass, is 



l{Mt,r)r) = ^Jlik)dk. (2.15) 

ko 

Let us now return to the discussion of the magnetic field energy and of the mag- 
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netic field structure, the topics which are of most importance in the kinematic dynamo 
theory. The evolution of the magnetic field energy spectrum in a turbulent conduct- 



ing plasma was first independently studied by Kazantsev ||2^ , and by Kraichnan and 
Nagarajan |^ in 1967. Kraichnan and Nagarajan were interested in the limit of 
small Prandtl numbers, which is not the case in galaxies and protogalaxies. Kazant- 
sev obtained the correct equation for the evolution of the magnetic energy spectrum 
on small subviscous scales in the limit of large Prandtl numbers. He assumed that 
the turbulent velocities are 5-correlated in time, (U(t)U(t')) oc 5{t' — t). In 1982 
Vainshtein derived a universal equation for the longitudinal correlation function of 



magnetic field |5T|]. He also proved, making only very general assumptions about 
the turbulence, that the magnetic field is exponentially amplified by the turbulent 
motions (i. e. that the fast dynamo inductive action does exist). It is interesting that 
in spite of more general assumptions Vainshtein used, his equation for the magnetic 
energy spectrum on subviscous scales basically coincides with that of Kazantsev 0. 

The complete correct theory for the evolution of the magnetic field energy spec- 
trum and for the growth of the total magnetic energy in the limit of large Prandtl 
numbers was first developed by Kulsrud and Anderson in 1992 They used the 
kinematic dynamo model and solved equation (|2.1j ) by the quasilinear expansion in 
time (by iterating it twice in time). Below we list the main results obtained by Kul- 
srud and Anderson, because in this thesis we will use calculational methods similar 
to those that they used, and the results we will obtain reduce to their results if the 
plasma is not magnetized (when the magnetic field is so weak, that cjjtj <^ 1, and the 



kinematic dynamo model is valid, see Figure |1.5| ). 

The equation (25) of Vainshtein and the equation (31) of Kazantsev are the same in the 
absence of resistivity, except there are different constant factors in front of their right-hand-side 
parts. The equation (2.25) in this thesis, which was derived by Kulsrud and Anderson would 
also be the same if converted to the differential equation for function k~^M . 
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Kulsrud and Anderson found that the ensemble averaged total magnetic energy 
per unit mass, E = (i?^)/87rp (p is the plasma density), grows exponentially in time: 

f - 2,.., (2.16) 
7o = \j:k'Jok. (2.17) 

They defined the magnetic energy spectrum M{t, k) as 

M{t,k) 4^(^)7 e {\B{t,k)\'} d% (2.18) 



Airp V27r, 

where the integration is carried out over all directions of k = k/k. Here, we consider 
that the function Bk(t) is the Fourier coefficient of the magnetic field B, while B(t, k) 
is the appropriately defined function continuous in k, so that 

r 3 °° 

5„(t,r) = E^k«(t)e*^^= (— ) I Mt,k)e''^''d'k (2.19) 



^ —CO 



(see definition of the Fourier transformations and coefficients in Appendix 0). The 
averaged total magnetic energy per unit mass is obviously equal to 



1 ~ 



^ = - / M{t,k)dk. (2.20) 







Kulsrud and Anderson derived the equation for the evolution of the magnetic 
energy spectrum M{t, k) 0, They called it the mode coupling equation^ which, in 



the absence of helicity, is 
dM 



dt 



I KJk, k')M(t, k') dk' -2 — k'^MU, k), (2.21) 
J An 
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where Q 



k" 



A I ^ \ I , k"^ + k''^ -kk' cose . , , 



2rf ^1)7^0.. 



(A;2 + A;'2 - 2kk' cos i 



A/2 



(2.22) 
(2.23) 



and 6' is the angle between k and k', see Figure [4.1| . The constant %o is the turbulent 
magnetic diffusivity: 



k 



1 



(2.24) 



Kulsrud and Anderson found that the magnetic energy cascades down to very small 



subviscous scales via the mode coupling equation (|2.21| ). In particular, this mode 
coupling equation greatly simplifies for scales much less than the viscous cutoff scale, 
k ^ ky^ and becomes a simple differential equation, homogeneous in /c. 



(2.25) 



If M(t, fcref) is known as a function of time at some reference wave number k = k^-et, 
then the solution of (|2.25|) is 



M(t, k) 



M(t',M Goik/h,i,t-t') dt' 



(2.26) 



where the Green's function Go{k,t) is 



5 \^/^ ^^^^Infc ^(3/4)^^j^_5ij^2^/4^^^ 



(2.27) 



^ Note, that Kulsrud and Anderson used a slightly different definition of the continuous 
Fourier transformation in time. As a result, they have additional factors 27r in their formulas. 
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A "signal" M{t, /cref) at zero time will increase exponentially as e^^^'^^"'°'^ and will extend 
down to the scale kpeak ~ e'^°*/cref, where /cpeak is the peak of function kGo{k,t). As 
a result, the magnetic energy tends to quickly propagate to very small subviscous 
scales d, 13 ! 

The next important question concerns the magnetic field geometrical structure, 
which is produced by the kinematic dynamo action. This structure has recently been 
extensively studied. It was first found that the mean square curvature of the field 
grows exponentially in time, oc e^^^/^^'^°*, and quickly becomes much larger than the 



square of the viscous wave number k^, |3T||. Even if the curvature is initially zero, 
it starts to grow because of the second order spatial derivatives of the turbulent 
velocities. Schekochihin et al. then found the actual probability distribution for the 
magnetic field curvature and showed that the curvature stays of the order of the 
viscous wave number in most of the volume, but becomes extremely large in a small 
fraction of the volume [^. Moreover, it turns out that there is an anti-correlation 
between the value of the curvature and the strength of the magnetic field! Namely, 
the magnetic field is weak in those (small) regions of the volume, where the field 
curvature is large, and the field is strong in those (large) regions of the volume, where 
the curvature is weak. Thus, the curvature stays finite ( ~ the smallest eddy size) 
over the bulk of the volume, while the magnetic field energy becomes concentrated 
on very small subviscous scales, as found by Kulsrud and Anderson. Therefore, the 
field has the folding structure with the wave number perpendicular to the field lines 
much larger than the wave number parallel to the field lines, k^ ^ kn ~ k^, as 



demonstrated in Figure This folding structure was first suggested by Cowley in 
1999 |TT]] on intuitive grounds and was later supported by MHD numerical simulations. 
The folding nature of magnetic fields is also theoretically supported by the Barnes 
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Figure 2.1: The folding structure of the magnetic fields produced by the kinematic 
dynamo (for simplicity shown in two-dimensions). The field is nearly straight and 
strong in Region I. The field is very curved but weak in Region II. 
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collisionless damping which states that small perturbations of the field (MHD 
waves) are strongly damped by collisionless energy transfer to the electrons The 
folding structure of the magnetic field lines is very important for our calculations of 
the magnetized dynamo! 



2.2 The Magnetized Turbulent Dynamo 

In a protogalaxy the magnetic field is quickly amplified by the kinematic dynamo 
action. When the field strength becomes of the order of ~ 10^^^ G or higher, the 
ion cyclotron time in the magnetic field becomes less than the ion collision time, 
ujiti ^ 1, (and correspondingly, the ion gyroradius becomes shorter than the ion 
mean free path). As a result, the plasma becomes strongly magnetized, and the 
dynamics of the turbulent velocities on the viscous scales starts to be controlled 
by the magnetic field. These are the scales on which the velocities most rapidly 



amplify the field (see footnote on page 11^). As a result, the main assumption of 



the kinematic dynamo theory, the absence of the back reaction infiuence of the field 
on the turbulence, becomes invalid. The turbulent dynamo becomes what we call the 



magnetized turbulent dynamo, see Figure |L5 



It is important to know that the magnetic Lorentz forces can be neglected as long 
as the magnetic energy stays less than the kinetic energy of the smallest turbulent 
eddies, i. e. as long as the field strength stays less than ~ 10^ G in a protogalaxy. 
During the magnetized dynamo stage in the protogalaxy, when tUjtj ^ 1, but the mag- 
netic energy is small, see Figure |L^, the magnetic field strongly affects the turbulent 



^ It works as follows. Sinusoidal MHD waves slightly perturb the field strength. As a result, the 
perpendicular kinetic energy of electrons oscillates in time. However, because the electrons refiect 
from magnetic mirrors, a part of the fluid kinetic energy is irreversibly transfered into the parallel 
kinetic energy of electrons (it is the second order effect). Thus, the electrons are heated, and the 
waves are damped. 
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motions on the viscous scales by changing the viscous forces, while the Lorentz forces 
are negligible! Indeed, under the condition that Uiti ^ 1, the transport of both the 
transverse and the longitudinal components of ion momentum in the direction per- 
pendicular to the magnetic field lines is inhibited. The transport of the transverse 
component of ion momentum along magnetic field lines is also inhibited, because gy- 
rating ions quickly lose their "memory" of the transverse ordered velocity in a time 
equal to the gyrorotation time 0. The transport of the longitudinal component of 
ion momentum along the magnetic field lines is the same as in the absence of the field. 
As a result, the viscous forces acting on turbulent velocities V in an incompressible 
fully ionized plasma are determined by the Braginskii viscosity stress tensor |^ 

vra/3 = -i^{3babp - Sa^)bf,KV^^u , (2.28) 

where h = B/B is the unit vector along the magnetic field. Note, that this stress 
tensor depends on the field unit vector b, but is independent of the magnetic strength 
i? = |B|, i. e. the viscous forces depend only on the magnetic field direction (as long 
as the plasma is strongly magnetized, and Uiti ^1). 

The MHD equations for the turbulent velocities V in an incompressible strongly 
magnetized plasma are |29| 



= - + fa + MbahbAy,^,^'),/^ - iyaVi3),fS , (2.29) 

Va,a = 0, (2.30) 

where f is the force driving the turbulence. In particular, it can be thought as a 
force, driving a turbulent eddy, that comes from larger eddies. The hydrodynamic 
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pressure P' can be determined by the incompressibility condition ( p.3(]| ). To obtain 
the second hne in equation ( |2.29| ), we use formula ( p.28| ) for the viscous stress, and 
we incorporate the isotropic part of the viscous stress into the pressure P" (for the 
same reason the force f can be assumed to be solenoidal). 

It is difficult to solve equations (|2.29|) and ( p. 301) directly because they are non- 
linear and they include the complicated Braginskii viscous forces. We also do not 
know the exact expression of the driving force f , but its statistics is the same as for an 
unmagnetized plasma. Therefore, let us proceed as follows. We introduce subsidiary 
incompressible turbulent velocities U which, by definition, satisfy equations 



dtU^ = -p;;: + u + ^uAu^-{UaUp),p, (2.31) 

Ua,a = (2.32) 



{AUa = Ua,/3i3 , and the pressure P'" can be determined by the incompressibility 
condition). Let us analyze and compare equations (|2.29|) and ( p.31| ). 



First, let us for the moment formally average the Braginskii viscosity term 
^v{hj) pb JjyV^^y) in equation ( p.29| ) over all directions of an isotropic magnetic 
field, i. e. let us apply formula {hj)phj)^)^ = (l/15)(5„/35^i, + S^^^Sp^ + 5au5p^) 
to it. Then this term reduces to (1/5)z/AVq, which coincides with the isotropic 
viscosity term in equation ( |2.31| ). In other words, z/efr = (1/5)^ could be consid- 



ered as an effective reduced viscosity for an incompressible fully ionized plasma 
in the presence of a magnetic field that is isotropically tangled on subviscous 
scales. 



Second, note that equations ( p.29|) and (|2.31| ) have the same driving force f . 



By taking the driving force to be the same, we assume that this force comes 
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from larger turbulent eddies. These larger eddies are on scales larger than the 
viscous scales, and therefore, these eddies "do not know" whether the viscous 
forces are of the Braginskii type or of the standard isotropic type. 



Third, note that equation ( p. 311) is a familiar hydrodynamic equation with a 



standard isotropic viscosity term. However, is has a reduced molecular vis- 
cosity, (l/5)z/ instead of u. Therefore, we can suppose that the solution of 
equation ( |2.31|) is the Kolmogorov turbulence with the effective reduced viscos- 
ity 



Z/eff = -V. (2.33) 



As a result, we suppose that the statistical properties and the kinetic energy 
spectrum of the subsidiary velocities U are given by formulas ( p^ , ( pTSf ), ( p.lO|) - 



( p.l5|) , where the Reynolds number R and the viscous cutoff wave number A;„ 



are now determined by the effective viscosity, z/eg = (l/5)z/, i. e. i? ~ VoL/vcE = 
WoL/p, K ~ R^^'^ko = {5VoL/uf/^ko. 

Fourth, we can say that because of the presence of the magnetic field, the Bra- 
ginskii viscous forces "are doing worse" at dissipating the turbulent motions, 
as compared with the standard isotropic viscous forces in a field-free plasma. 
This is why the effective viscosity in equation ( p^.31D is reduced by a factor 



1/5. Thus, we can draw the important conclusion that the spectrum of tur- 
bulent velocities in a strongly magnetized plasma extends to a smaller scale, 
27ik~^ ~ (5Vo-L/z/)^^/^L, as compared to the cutoff scale in an unmagnetized 
hydrodynamic plasma, 27ik~^ ~ (VoL/z/)~^/^L! [] 



^ This conclusion is important because it means that in a magnetized plasma the smallest tur- 
bulent eddies, which are located on the viscous scales, have turnover times shorter than those of the 
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Now let us subtract equations ( p.31| ) and ( |2.32| ) from equations ( |2.29| ) and ( p.30| ) 
to eliminate the unknown driving force f, and let us introduce the hack-reaction 



velocity V =^ V — U. We have 



dtV 



- {VaUp + U^Vp + v^Vp)^p , (2.34) 
0, (2.35) 



where the pressure P = P" — P'" . Here and below we use the following symmetric 
tensors (which vary in time and in space): 



def 



babfsb^bs , 



— babijbj , 
baf3 = bj)^ 



(2.36) 
(2.37) 
(2.38) 



Velocity V, which satisfies equations ( p.34|) and ( p.35|) , can be considered as the 
correction to the Kolmogorov velocity U. This correction results from the strong 
influence of the magnetic field direction la on the turbulent motions through the 
Braginskii viscosity tensor (|2.28| ). The evolution of the magnetic field B is given by 



the following standard MHD equation 



(2.39) 



where the plasma velocities \ = \J + V are incompressible [compare with equa- 

smallest eddies in an unmagnetized plasma. As a result, the magnetic field amplification rate in the 
magnetized plasma is expected to be higher, compared to that in the unmagnetized plasma. See 
also footnote on page fsl. 
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tion and we neglect resistivity. Consequently, the equations for the magnetic 

field squared, B^, and for the magnetic unit vector h are 



(2.40) 



dtb, 




(2.41) 



In the end of this section, it should be noted that the magnetized turbulent dynamo 
does exist in protogalaxies, but does not in galaxies. The reason for this is that, in 
contrast to temperatures in protogalaxies, temperatures in galaxies are low enough 
(T < 10^ K) for a considerable number of neutral particles to be present. These neutral 
particles have mean free path much longer than the ions have. Therefore, the viscosity 
in galaxies is dominated by the neutral particles, and the viscous forces are given by 
the standard familiar formula, z/nAV (here is the neutral viscosity). As a result, 
in galaxies the turbulence is Kolmogorov and the kinematic dynamo theory can be 
used 0. 



* Temperatures can be very high in the coronal regions of galaxies, and in these regions the 
magnetized turbulent dynamo can be present. 
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Chapter 3 



Statistics of Turbulent Velocities 



in Strongly Magnetized Plasmas 



The statistics of the turbulent velocities in an unmagnetized plasma is given by for- 
mulas ( p.4|) and ( p.5|) . In this chapter we derive similar equations for the statistics 



of the turbulent velocities in a strongly magnetized plasma by making use of the 
quasi-linear (up to the second order) expansion in time, to solve the MHD equations. 
The final results of the derivations are given by equations (|3.8CI|) - (|3.83|) . We then use 



these results in Chapter ^ to derive formulas for the evolution of the magnetic field 
total energy and energy spectrum in the magnetized turbulent dynamo theory. 

3.1 The Quasi-linear Expansion of 
the MHD Equations 

Let us assume that we know the magnetic field at zero time, B|^^q = °B(r) and 
b|j^Q = °lo(r), and that the back-reaction velocity V is initially zero, I^Ij^q = ^v{r) = 
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0. ^ We advance the magnetic field and the back-reaction velocity to some future time 
t > by the nonlinear terms, i. e. by integrating equations ( |2.34| ), ( p.39| ) and ( p.41| ) 



twice in time. This quasi-linear expansion procedure is similar to the calculations of 
Kulsrud and Anderson Considering t as the expansion parameter 0, up to the 



second order, we have 



B(t,r) 


= °B(r) + iB(t,r) + 2B(t,r), 




b(t,r) 


= °b(r) + ib(t,r)+2b(t,r). 




V{t,v) 








V(t,r) 


= ^V(t,r) + 2V(t,r) = 







(3.1) 
(3.2) 
(3.3) 
(3.4) 



where V = U + 1^ is the total fluid velocity, and the incompressible Kolmogorov 
turbulent velocities U are considered to be given and to be of the first order 
All first and second order quantities, except U, are initially zero, e. g. ■'^B(0, r) = 
and ^V{0, r) = 0. 

Now, we substitute the above expansion formulas into equations ( |2.39|) , ( |2.41|) and 



into equations ( p.34| )-( p.37|) . We find that the zero order equations are 



dt'^B^ = 0, (3.5) 

dt% = 0, (3.6) 

''Va = 0, (3.7) 

V, = 0, (3.8) 

%al3y = %%%, (3.9) 



^ A nonzero initial back-reaction velocity would lead to transients which would be dissipated 
anyway. 

^ If we would like to be more formal, then we need to introduce a dimensionless variable ^ = t/ At, 
and to consider Ua^f^At and Va^f^At, which are dimensionless, as the expansion parameters. 
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'^b^p.u = %%%%; (3.10) 
the first order equations are 

dt = ^Va,l3^Bi3 — j3^B ^,13 , (3-11) 

dt^h^ = V„ A - V;3,^X/37 - , (3.12) 

dt^Va. = -'P,a + 3z/(X/3^.'l^^,.),/3 + 3z/(X^^,t/^,,),^-iz/At/a , (3.13) 

o 

'Va,a = 0, (3.14) 

V„ = U^ + 'Va, (3.15) 

^b^p-y = + + , (3.16) 

= + + + (3.17) 

and the second order equations are 

dt ^Ba = ^Va,l3^Bp + '^Va,/3^Bi3 — ^Vp^Ba,/3 — "^V p^B , (3.18) 

- 'Vf,%,fs-^Vp%a,f3, (3.19) 

- 'v^Up-Ujvp-'vJvp)^^, (3.20) 
X„ = 0, (3.21) 

V„ = V. (3.22) 



Here, of course, the pressure is also expanded, P = + ^P, and it is completely 
determined by the incompressibility conditions ( p.l4| ) and ( p.21| ). 
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3.2 The Time Laplace Transform Solution 
for the Velocities 



Let us first solve the first order equations ( |3.13| )-( |3.15| ) to find the first order back- 
reaction velocity and the first order total velocity "''V. It is convenient to Fourier 
transform these equations in space, r — > k, by making use of the discrete Fourier 
transformation, and to Laplace transform them in time, t — > s, by making use of the 
continuous Laplace transformation (see Appendix |^) . We have 



I3fj,u 



k' 

k"= k-k' 



+ 3uikf3 Yl ^KUk'fMisfbl,"a/3^,u + \j^k'^Uka{s), (3.23) 
k' ^ 



kjvi^ais) = 0, (3.24) 

'V^a{s) = U^a{s) + 'v^ais), (3.25) 



where ^V^a, ^Pk, ^"t^kv, Uy^a and ^V^a are the Fourier-Laplace coefficients, and 
^b\^"oti3^iv is the Fourier coefficient i^hi^nap^u is constant in time because it is of the 
zero order). In the derivation of equation ( |3.23| ) we also make use of the zero ini- 



tial condition, '^v\^^Q = [see equation ( [A.ll| ) in Appendix Now, we multiply 



equation ( ^.23| ) on the left by tensor = 6^a — k^k^, which is perpendicular to 
the unit vector k = k/A;, see equation ( p. 61) . The pressure term goes away. Using 
the incompressibility condition ( 3.24j) , interchanging indices, and using the symmetry 



property of tensor °&k"a/3/i!/ over indices a, /?, fi, and u, we obtain 

°-^ka,k'/3"'^'^'k'/3 = ^^ka , (3.26) 

^Mkayp = s(5k„,k/^ + 3z/(5^^fc^C^k"7M^'^ , k" = k - k', (3.27) 
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ko 



-°-Mk«,k'/3 + (^S + 5ka,k'/3 t/k'/3 . (3.28) 



Here we use convenient matrix notation, so that summation is implicitly assumed 
not only over repeated spatial indices but also over repeated wave numbers [e. g. a 
summation over (3 = x,y,z and an infinite summation over all discrete values of k' 



are assumed on the left-hand-side of equation ( |3.26| )]. Function 5ka,k'/3 = ^k,k"^a/3 is 
the Kronecker tensor, it is equal to unity if a = /3, k = k', and is zero otherwise. The 
matrix operator 'lMka,k'/3(s) is of the zero order, while "the driving force" \7-'ka(s) is 
of the first order. 

If there exist matrix 'iMkQ.k'/? ^^^X is inverse to the matrix 'iMka.k'/?, then the 
solution of equation (|3.26|) is 



^v^a = °-Mkik'/-^k'/3 . (3.29) 



Using this formula and equations ( p.25| ), ( |3.28D , we obtain the Fourier-Laplace coef- 
ficient of the first order total velocity -"-V, 



'^ka(s) = f/ko + 'C'ka ={s + °-Mkik'/3^k'/3(s) . (3.30) 



Now, let us consider an important case for which we can invert matrix ^A^ka.k'/?, 
given by formula ( p.27| ). Let (°B ■ V)°B = 0, i. e. let the magnetic field at zero time. 



°B, vary only in the direction perpendicular to itself. This is equivalent to the case 
when the magnetic field lines are initially straight, the initial curvature of the field is 
zero, and ^hap = ^hj^h^ = const. This model of straight magnetic field lines is not as 
artificial as it seems at first glance, because of the following arguments. Remember, 
that in the kinematic dynamo case, the magnetic field lines have the folding pattern. 
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shown in Figure |2.1| . This folding pattern imphes that in the bulk of the volume the 
field is strong and has small curvature, k± ^ ~ k^, while in a small fraction of the 
volume the field is weak and curved, k± ~ k\\ ^ k^. The regions of weak and curved 



field. Region II in Figure |2.1| , can be disregarded as long as we consider the volume 



averaged magnetic field energy spectrum and are not interested in the field curvature. 



As for the regions of strong magnetic field with small curvature. Region I in Figure ^ 
the field in these regions can be well approximated by our straight magnetic field lines 
model on scales k which satisfy k > k\\ kj^. However, it should be made clear that 
we make an assumption that in the magnetized turbulent dynamo case the magnetic 
field has a folding structure similar to that in the kinematic turbulent dynamo case. 
This assumption is our first working hypothesis. We do not prove it, but there exist a 
simple argument by contradiction, supporting it, which is as follows. If the magnetic 
field were not folded, i. e. if k±_ ~ k\\ > ki, everywhere in space, then the field would 
be isotropically tangled on viscous and subviscous scales. In this case the Braginskii 
viscosity could be averaged over the field directions, and it would probably rapidly 
reduce to the isotropically averaged viscosity u^q = z//5, see equation ( p.33|) . As a 
result, in this case the magnetized dynamo would probably possess the properties 
of the kinematic dynamo with this effective viscosity, and would develop the folding 
structure of the magnetic field. 

In the case when ^bap = const, we have °6k"7/i!//3 = Si^" fi^b^'^b^'^bu'^bp, and ma- 
trix ( p.27| ) reduces to 

%fka,k'/3 = 4,k' [sSaf, + Su^fl^S^fb^S] , (3.31) 

II %-k. (3.32) 

Here fx is the cosine of the angle between vectors "^b and k. The matrix '!Mka,k'/3 is 
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now diagonal in k, so we can easily invert it: 



''^'kQ,k'/3 



Jk,k' 



_ 3z/fcV(l -^^) 



s + 3z/A;2/i2(l - /i2) 1 - /i2 



(3.33) 



Next, we substitute this inverse matrix into equation ( p.30| ) and obtain the Fourier- 
Laplace coefficient of the ffist order total velocity 



ka^ 

s + Vt 



kal 



a7 



s + 2n 1 - /i2 



f/k/3(s) 



(3.34) 
(3.35) 

(3.36) 



where, we introduce the following notations for the viscous damping frequencies 



Q =^ -uk 



2 7 2 



def 3 



15^ 



-z.A;V(l-/i') = yf^/^'(l-/i') 



(3.37) 
(3.38) 



The frequency Cl represents the effective averaged rate of the Braginskii viscous dis- 
sipation, see equation ( |2.33| ). The frequency Q depends on /i2 = (°bk)2, and this 
dependence reflects the anisotropy of the Braginskii viscous stress tensor. Note, that 
fl is equal to Q averaged over fi. 

Now, we calculate (^Vkait)) and {^VkaityVk'^it')), which are the ensemble av- 
erages of the V's over all possible realizations of the turbulent motions. For this, 
ffist, we need the ensemble averages of the U's. We use the Laplace and Fourier 
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transformation formulas (see Appendix^), and equations (p.4|), ( p.5| ) to obtain 





oo 





oo oo 



{U^aioo)) e-'"^' duj = 



(3.39) 



{UUt)Ui.'f3{t'))e-''-''''dt'dt 

] oo oo 

e-^'-^'Ut'dt ^ f f {UMU^'pi^')) 







27r 

— oo— oo 
oo oo oo 



^ I -Lkdu I I e~''~'^'e-''+'^'dt'dt 



1 



-oo 

oo 



duj 



2tx J '^^ (s + iuj){s' — iuj) 



(3.40) 



Second, we use these formulas to calculate the ensemble averages of the two modes 
of the velocity Fourier-Laplace coefficient in equation (|3.34|) , ^V]^^{s) and ^V]^^{s). 
Using equations ( p.35| ) and ( |3.36| ), we obtain 



''v'Usfv'i',{s)) 



0, 
0, 
0, 

5k',-k 



ss' 



1 - 



X 



Jivk duJ 



2n J {s + iiu){s' — iuj) 



X 



Juk du 



2tt J (s + iu)(s' — iuj) 

-oo 



(3.41) 
(3.42) 
(3.43) 



(3.44) 



(3.45) 
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Third, using these last equations and equation ( p.34| ), we easily obtain 



(3.46) 



~JI , , 1 ~ /' 



where 



5k',-k — y JujkduJ 

— oo 



Hl{s;uj,0)Hl{s';-uj,0) 6^0- 



+ Hl{s]uj,VL)Hl{s']-u,VL) 



xl. x± 



(3.47) 



P + n 



{p + iqi){p + 2q2] 



(3.48) 



The inverse Laplace transformation p — t of function Hl{p] qi, ^2) is 



c+ioo 



HLit,qi,q2) 



2m 



Hl{p] qi,q2)e^^dp 



c—tca 
c+ioo 



2m J {p + iqi){p + 2q2) 

c— «oo 

{n - iqi)e-"^'^ -{Q- 2^2)6-^92* 
2q2 - iqi 



(3.49) 



where we carry out the integral by closing the integration contour in the complex 
plane and by calculating the residues at p = —iqi, —2q2. Next, we apply the inverse 
Laplace transformations s t and s' t' to equations ( |3.46| ) and ( p.47[ ). We finally 
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obtain the desired results 



where function Hiij:, qi, ^2) is given by equation ( p.49| ). As one might expect, the en- 
semble averaged first order turbulent velocity, {^Vi^a{t)), is zero 0, while the averaged 
velocity correlation tensor, (^V^ka(^)^^k'/3(^')) is not zero because it is of the second 
order. 

So far we have considered only the first order velocity, -"-V, and have found its 
statistics, given by equations ( |3.50| ) and ( ^.511 ). In order to find the second order 
velocity, ^V, we need to solve the complicated second order equations ( p.20|) -( p.22| ). 
Fortunately, we will need only the ensemble average of the second order velocity, 
(^Vkait)). It turns out that in our case of a straight initial field, °6a/3 = const, which 
we consider here, this average is zero, 

CVUt)) = 0, (3.52) 

see Appendix ^. The reason for this simple result is that ^6^/3 is constant in space, and 

This is because there is no preferred direction. Of course, there is a preferred axis in space, 
which is along the magnetic field unit vector. However, there is no preferred direction because the 
Braginskii viscous stress tensor ( [2.2^ ) is invariant with respect to a substitution b — > — b (gyrating 
ions "do not care" about the exact direction of b). 



(3.50) 



5k',-k — / Ju)k du 



x± x± Oi Oi 

H,{t; 0)H,{t'- -u, 0) I 6^, - '"Y: ^2 

x± x± oi Oi 

+ H,{t; uj, m^it'; -co, fi) "^^^ I ^ 

1 — /i^ 



(3.51) 
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the hydrodynamic turbulence U is homogeneous. Therefore, the ensemble average of 
the term in the parentheses in equation ( p.20| ) is constant in space, and its spatial 
derivative is zero. As a result, the ensemble averaged velocity (^Vkait)) = (^Vkait)) 
is zero (see Appendix for the proof). 

Now, let us return to the first order velocities, -"-V. According to equation ( p.50|) 



the ensemble average of """V is zero, ("""V) = 0. Let us calculate the ensemble average 
of ■'■V squared. We have 



(V) = (V„(t,r)V„(t,r)) = 5:(irk„(t)^V^k'„(t))e^('^+''> 



k,k' 

T.cyutYv-ut)) = E(rvk(t)r) . (3.53) 



Here we use the fact that C'Vka^Vk'is) = (^^ka^'^V^-k'/s) ^ ^k',-k, see equation (|3.51| ). 
Using equations ( p.49|) and (|3.51| ), we obtain 



00 

(rVk(t)P) = ^ f J^kduj[HL{t;u,0)HL{t;-LO,0) + HL{t;iJ,Q)HL{t;-LO,Q)] 
Ztt J 



—00 



00 , 
1 r ^ , ^9 1 — cosu;t - smut 
— / J^kduj\2Q^ + 2Q + 1 



-00 



[Q - (Q - 2Q)e-^^f 2o^l-coscjt 
+ 2(n-2n)e--4^ + ^^|. (3.54) 



Let us analyze this last equation. The width of function J^k in ^ space is about 
1 1 T{k), where r(fc) is the decorrelation time of hydrodynamic eddies on scale 27[k^^, 
e. g. see equation ( p.lO|) . We have an obvious estimate l/r(A;) < 1 /t(A;,^) ~ i^csk^, 



where ki, is the viscous cutoff wave number, and is the effective viscosity (|2.33| ). 



In order to make an estimate of the right-hand-side of equation (|3.54| ) , let us for the 



46 



moment consider the case when the angle between k and b is neither very close to 
zero or vr, nor very close to ±7r/2, so that 2Q = 15Cl^'^{l— fi'^) is not small (/i = k-°b). 
Then, in the limit t ^ r we have Qt ^ 1, Clt ^ 1 and uj^^ sin cut ~ n6{uj). Therefore, 
we can drop the exponential terms in equation ( p. 54]) , and the first term in the 



brackets {...} in this equation gives the main contribution to the integral over u. As 



a result, equation (|3.54|) reduces to 

(rVk(t)H ^ Jo^nH ^ t, t:$>T. (3.55) 

The velocity modes, which have the angle between k and °b not close to 0, vr, or ±7r/2, 
grow in time 0. Because all terms under the sum in equation (|3.53|) are nonnegative. 



the ensemble average of the first order velocity squared, ("""V^), also grows in time. 
It is clear that in reality this growth can not be as unrestricted, as it appears to 
be, according to equations (|3.53| ) and (|3.55|) . In the next section we will discuss this 



growth in more details. We will see that this growth is a feature of the anisotropy 
of the Braginskii viscous forces and of our quasi-linear expansion procedure, and is 
indeed restricted. 

3.3 The Effective Rotational Damping 
of Velocities 

Let us try to understand the reasons of the growth of the ensemble averaged first 
order velocity squared, ("""V^), which we found in the previous section. First, refer 
to equations ( p.43| )-( p.45| ) and equations ( p.47| ), (|3.51| ). From calculations of the first 



^ If the angle between k and °b is 0, tt, or ±7r/2, then = 0. In this case equation ( 3.54 ) 
obviously reduces to (|''"Vk(t)P) = 2Jofcf2^i oc t, and the growth is even faster. 
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order velocity growth, which follow equation ( p.51| ), it is clear that in the case when 
the angle between unit vectors k and °h is not equal to 0, vr or ±7r/2, the growth of 
■'■V(t) = "'"V'(t) + ^\"{t) happens due to the growth of the ^V^ka(^) modes, while the 
^V]f_ait) modes do not grow unrestrictively 0. Next, let us refer to equations ( p.35D 
and dpI) for ^V^^ and ^V^^. On one hand, we have 



kJKais) = 0, (3.56) 
L'Kais) = 0, (3.57) 



as it should be because plasma velocities are incompressible. On the other hand, we 
have 

= 0, (3.58) 

= ^^'^b^UUs). (3.59) 

Thus, the growing velocity modes -"-Vj^ are perpendicular to both vectors k and °b (of 
course, all modes must be perpendicular to k because of the fluid incompressibility 
condition). At the same time, the other, non-growing modes, """Vj^, have nonzero com- 
ponents along the initial magnetic field unit vector °b. In addition, equation ( p. 431) 
means that V,^ and V,^ are perpendicular on average. 

These results can be understood by an independent physical argument that ap- 
peals to the Braginskii viscous dissipation process. The viscous dissipation into heat 



^ This can also be seen because of the following simple arguments. The poles of function 
(^^ka('5)^^k'/3(*')) ^^'^ s — 0, s' = 0, s = —iuj and s' = iuj, see equation (3.44). All these poles have 
nonnegative real parts, and there is no viscous damping of (|"''Vi(.(t)p). On the other hand, the poles 



of function ("^V^ka('S)^^k'/3(*')) ^"^^ * = — 2f2, s' = — 2f2, s = —iuj and s' = ioj, see equation ( 3.45 ) 
Two poles are negative, and (|''"Vi(.(t)p) is viscously damped with a rate ~ 51. 
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IS 



Qvis = T^apV^fi = - 3z/ {Wa/h3pf = - 3z/ ('V''%3pf , (3.60) 



where we keep only the first order terms in the parentheses, and make use of equa- 
tion ( |3.58| ). Using this formula, the inverse Laplace transform of equation ( p.45[ ) 



[which is the second term in the square brackets in equation ( p.51| )] and defini- 
tion ( p.38| ), we obtain the ensemble averaged dissipation 



k,k' 

k 

1 °° 

= -Y.'^^T^ [ J^kHL{t;u;,n)HL{t';-uj,n)duj . (3.61) 
iT 

" — oo 

We see that first, only ■'■V,^ velocity modes are dissipated by the Braginskii viscous 
forces, and second, the viscous dissipation is proportional to 2Q = 3z/A;^/i^(l — /x^). 
As a result, we can conclude that the growth of the first order velocity happens 
because the """Vj^ velocity modes, which are perpendicular to °b, are not damped by 
the Braginskii viscous dissipation process. Of course, in the degenerate cases, when 
k ± °b (i. e. = 0) or k || °b (i. e. 1 — /i^ = 0), both velocity modes ■'■vj^ and ^v'^ 
are undamped. 

It is clear that the growing velocity modes can not grow unrestrictedly. What are 
the mechanisms which stop this growth? To answer this question, let us again note 
that the only growing velocity modes are those which are perpendicular to b. All 
other modes are damped by the Braginskii viscosity. For example, in our quasi-linear 
expansion the zero order magnetic field unit vector °6 is assumed to be stationary. 
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see equation ( p. 61) , and therefore, the first order velocity mode "^V,^(t), which is per- 
pendicular to °b, can unrestrictedly grow in time. However, in reality, any velocity 
vector continuously rotates relative to the magnetic field unit vector. As a result, the 
growing modes do not stay perpendicular to h, and therefore, they do not grow in 
time forever. We can understand this rotation from different standpoints. First, in 
the laboratory reference frame the velocity vector and the magnetic field unit vector 
both rotate, but they must rotate differently because different "forces" are acting on 
them [in other words, because they satisfy different differential equations]. Second, 
we can go to a reference frame that has the origin at a given point of space and rotates 
together with the field unit vector at this point. In this rotating reference frame there 
exist Coriolis forces which act on velocity vectors and force them to rotate relative 
to the non-rotating field vector. These Coriolis forces are caused by the non-linear 
coupling of velocity modes via the inertial term (VV)V of the MHD equation (|2.29| ). 
Consequently, the rotation is produced by the non-linear coupling of velocity modes. 

Thus, a growing velocity mode, which is perpendicular to the magnetic field unit 
vector, rotates out of its initial direction. After this, the velocity mode is viscously 
damped, and stops growing. We call this process the "effective rotational damping" 
of velocity modes. Of course, the effective rotational damping operates only on the 
viscous scale, on which the Braginskii viscous dissipation is significant. On larger 
scales the viscous dissipation is small and the rotation of velocities does not make 
any difference. Now, let us estimate the effective damping rate of the growing velocity 
modes, which is associated with the rotational damping process. First, the square of 
the angular velocity of the rotation can be estimated as 

u^L-l^, (3.62) 
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where is the velocity decorrelation time on the viscous scale. The factor 1/3 in this 
equation comes from the fact that only one of the three angular velocity components 
contributes to the deviation of the velocity V from its original direction perpendicular 
to the magnetic field unit vector h. Namely, only the component along the vector 
product b X V contributes. Second, the typical Braginskii viscous damping rate 
is about 1/r,^ ~ T-'eak'l because of the definition of the viscous wave number ky, as 
the wave number at which the viscous dissipation becomes important. Thus, the 
typical angular velocity of the rotation of the velocity modes is less than the typical 
viscous damping rate. Therefore, we assume that after a time interval At the growing 
velocity mode rotates by an angle A0 relative to its original direction perpendicular 
to b, only the projection of the velocity mode on an undamped direction (which is in 
the plane perpendicular to b) survives, and all other components of the velocity mode 
are immediately viscously dissipated. As a result, the effective rotational damping 
rate fird can be estimated as 



dV IW V(cosA0-l) V V , ,2^t 

-dE ^ AT A^ 2At ^ - 2At - 

~ - ~ ^^efffc^V^^^ ~ h^^'^y' = - ^rd^'' (3-63) 

^^rd = ^i^eskl = j^ukl. (3.64) 



To obtain the last result in the first line of equation (|3.63|) , we use the random- walk 

approximation for the estimate of A0^. f\ We also use formula (|3.62|) for tOj-ot- 

As we said, the effective rotational damping is associated with non-linear coupling 

of velocity modes. On one hand, this non-linearity is hard to deal with directly. On 
° Basically, here we assume that 

^rot^u ^ -^/^ 1" this case we can choose such time 
interval At 3> that ~ tJiotAi ^1. As a result, on one hand, we can expand the cosine in 



equation ( 3.63| ). On the other hand, the turbulent eddies on the viscous scale are correlated on time 



intervals ^ ^ At, and we can use the random-walk approximation to estimate A(/)^ 
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the other hand, the rotational damping is physical and is very important, it has 
to be included into the theory in order to restrict the growth of the velocity modes 
undamped by the Braginskii viscous forces. Therefore, we are going to incorporate the 
rotational damping into our equations in a simple way as follows. First, we go back to 
the original equations for velocities, equations ( p.29|) and (p.34|) . We recall that the 
total velocity V is the sum of the Kolmogorov hydrodynamic turbulent velocity U, 
which is assumed to be given, and of the back- reaction velocity V, which we obtain 
by solving equation ( |2.34] ). The Kolmogorov turbulence is assumed to be stationary 
(and homogeneous). Thus, the growth of the undamped modes of the total velocity 

V = U + V happens because the corresponding modes of the back-reaction velocity 

V are undamped and grow [| The rotational damping restricts the growth of these 
back-reaction velocity modes on the viscous scale. As for the scales larger than the 
viscous scales, the turbulence is Kolmogorov on these large scales, V = U, and the 
back reaction velocities are zero [see the discussion in Section |2.2| , also note that the 
driving force ( |3.28D , which is oc /c^, becomes small on large scales]. As a result, on all 



scales we add the rotational non-linear damping into our equations for the velocities 
V. We do this by replacing the time derivative d/dthj d/dt + VL^^ in equation (|2.34| ): 



- {VaUf3 + UaVp + VaVp),fS , (3.65) 

v^^o^ = 0. (3.66) 

Once again, the purpose of this replacement is to incorporate the rotational non-linear 
damping, which stops the unrestricted growth of the velocity modes undamped by 



^ One can obtain the Fourier-Laplace coefficient of the first order back-reaction velocity, ^Vk, 
from equatio n (|3.30 ), and then derive formulas for the growth of (■'■'L'^(i)), which are very similar 



to equations ( 3.5E )~( 3^ ) 
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the Braginskii viscosity, into our equations in a simple way, by including the linear 
damping term Qj-dVa- This is our second working hypothesis. A possible check of 
it can be intensive MHD numerical simulations, including the Braginskii viscosity. 
However, such simulations are rather complicated |, and they are beyond the scope 
of this thesis. Note, that in equation ( p.65|) , the effective rotational damping term 
fij-d'^a is isotropic, and therefore, it damps not only the growing modes of V (which are 
undamped by the Braginskii viscosity) but all V modes. This is not a serious problem 
though, because the rotational damping is smaller than the Braginskii damping by a 
factor ~ 1/6 [see equation ( |3.64|) ], and our results should be valid within a factor of 



order unity. 

3.4 The Time Fourier Transform Solution 
for the Velocities 



In this section we carry out calculations similar to those of Section p.2| , again assuming 
our first working hypothesis that the initial field can be taken to be straight, ^bajs = 
const. We also again use the quasi-linear expansion ( |3.1| )- (|37^ ). The main difference 
is that now, instead of using equation ( p.34| ) for the back-reaction velocity, we use 
equation ( p.65| ) with the effective rotational damping ( |3.64[ ) in it. In the end of this 



section we finally obtain convenient formulas for the statistics of turbulent velocities 
in a strongly magnetized plasma. 

The first and the second order equations (|3.13| ) and (|3.20|) for the back-reaction 



^ In particular because the Braginskii forces are anisotropic and there are undamped velocity 
modes, it is not possible to drop the inertial terms in the MHD equations (i. e. not possible to assume 



the viscosity dominated regime). For example, if we drop the inertial term sSap in equation (3.31) 
then the matrix °A^ka,k'/3 will not have an inverse matrix. 
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velocity become 



luAU^ , (3.67) 
5 



- ^v^Up-Ujvp-^vJvp) (3.68) 



and they now include the rotational damping terms in them. All other first and sec- 
ond order equations (pll)-(pl2|), (^)- ( ^39]) and (lOTD-dO^) stay unchanged. 



In Section |3]^, to solve the first and the second order equations, we use the Fourier 
and the Laplace transformations in space and in time respectively. In this section it 
is convenient to use the Fourier transformations both in space and in time 0. We use 
the discrete Fourier transformation r ^ k in space and the continuous Fourier trans- 
formation t — ^ in time (see Appendix Comparing equations ( p.l3| ) and (|3.67| ), 



and equations ( |A.9| ) and (|A.12|) , we easily see that the double Fourier transforma- 
tion (in space and in time) of equation ( p.67| ), which is the first order equation for 
the back-reaction velocity Fourier coefficient ^€'k,o('^), coincides with equation ( p.23| ) 



where the Laplace variable s is replaced by —iuj + l^rd- As a result, there is no need 
to rederive all equations of Section p.2| . In most cases we can obtain the new formulas 
just by making the replacement s — > —iuo + i^^d in the corresponding formulas of Sec- 
tion |3.2| . However, let us be more general, and in the calculations below assume that 

the effective rotational damping rate VL^.^ is a function of = |k| and of (b • k)^, |3 

^ Using the Fourier transformation in time is actually better than using the Laplace transforma- 
tion. The reason is that in the case of the Laplace transformation we have to assume the initial 



condition on the back-reaction velocity (in Section 3^ we assume ""■I^Ij^q), while the Fourier trans- 
formation "docs not care" about the initial condition. As a result, we do not have to worry about 
any velocity transients, which are eventually damped away by the Braginskii viscosity and by the 
rotational damping. 

We choose ilj-d to depend on the square of b • k because the Braginskii viscosity is invariant 
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e. g. Qrd = ^rd(fc, At^) up to the zero order, where fi is given by equation ( |3.32| ). We 
can always substitute our simple estimate ( |3.64| ) for Qj-d into our final formulas. 

Equations ( p.34D -( ^.36[ ) for the Fourier coefficient of the first order velocity now 
become 



—iuj + Vt + fifd 

—iuj + r^rd 



.UJ) 



1-/^2 



-iLu + n + ^ifhfhp ^ 



(3.69) 
(3.70) 

(3.71) 



In these formulas the effective rotational damping rate, fird, is added to the viscous 
damping frequencies ( |3.37| ) and ( p.38| ). Thus, the overall damping is produced by both 
the rotational damping rate and the viscous damping (of course, there is only one 
physical dissipation process present in the fluid — the viscous dissipation). Now, we 
use formulas ( ^.4] ), ( |2.5| ) to obtain averages similar to those given by equations ( ^.411 )- 
(g). We have 







(^)) 


= 








(^)) 


= 0, 








= 






'K'pi 




— Jujk 5k',-k 








X 





rd J 



rd 



!± S.± Ol^ 



CKaiooYKpi^')) = ^.5k',-k5(^' + ^) 



u;^ + {n + ad)' 
+ {2n + fird)' 



(3.72) 
(3.73) 
(3.74) 

(3.75) 



under reflection b ^ — b. 
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X 



l-/i2 



(3.76) 



Using these equations and equation ( |3.69| ), we obtain equations, which are similar to 
equations ( p.46| ) and ( p.47| ) of Section 



CV^aico)) = 0, (3.77) 



X 



1-/^2 



1 - 



where function Hf{uj] qi, 52) is 



+ ql ^ _^ - qi 



uj^ + qi 



(3.78) 



(3.79) 



Next, we apply the inverse Fourier transformations uj t and uj' t' to equa- 
tions (|3T7|) and ( pTSI) . We have 



(^V^k.(t)) 

('^ka(t)'V^k'/3(t')) 



0, 



(^V^k.(t)^V^*_u',(i')) 



(3.80) 



5k' 



k',-k 



_ + ad, ad) 5i 



1-^2 



r± f± oe 0£ 

+ Hpit -t'-n + ^d, 2fi + ad) ! ' 

\ — jj,-^ 



(3.81) 
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where function 



1 °° 



-co 
oo 



— / JujkHF{uj;qi,q2) COS [uj{t -t')]duj (3.82) 



depends only on the absolute value of the time difference t — t' because J^jk and 
Hp{(jj; qi, ^2) both are even functions of u. 

Equations ( |3.80|) - (|3.81|) are similar to equations ( p.50| )-( |3.51| ) of Section p?2| , 



and they give the statistics of the first order turbulent velocities, which we will use 
in the next chapter. As, for the second order turbulent velocity, equation ( p.52| ) stay 
the same, 



(3.83) 



because the proof of this formula, given in Appendix does not involve any trans- 
formation in time. 

In the end of this section let us substitute equation ( 2.10|) for J^^ and equa- 
tion ( p.79|) for Hf{uj; qi, 52) into equation ( p.82| ), and obtain the following formulas 



1 J 

—00 



1 + 



uj"^ + ql 



Jok _|_ Jok Ql — (I2 



2t 2 g2(l + rg2) 



(3.84) 



It is interesting, that we can obtain equations ( |3.80| )-(3.82) from equations ( 3.49| )-( ^.51[ ) by 
first, making replacements — > H + fird, 2(72 — >■ 2g2 + f^rd in equation ( 3.49| ), and second, by taking 
a limit t — *■ cxd in it. This is not surprising, because as time goes on, the velocity transients, which 
"remember" the initial conditions, are damped away by the Braginskii viscosity and by the rotational 
damping. 
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t t' 







t t 



J dt'J Hpit' - t")df 







1 7 T 
i / J( 



Ok 



2tt J 1 + r2cj2 

oo 



1 + 



u;2 + ql 



sinci;t 



duj 



Jok/^ ^-^tlT\ . Jok Ql — qi 



2 _ „2 

7 7 2 2 7 2 



1 - 6-92* 



-^/r^ 



2 ' 



1 dt' jHpit' -t")dt" = ^ j 



1 7 J. 



OA; 



TT J 1 + r^a;^ 

00 



1 + 



2 ^2^ 

?! - qI 

+ g2 



(3.85) 



sin^ (cut/2) 



t - r(l - e'*/") 



+ 



Jok Ql — Q2 



q2t -1 + e 

qI 



-q2t 



2 l-r2g2 

Jpk , Jok_ Ql — Q2 , _ Jpk Ql , 

2 2 g2 - 2 g2 ' 

t t' 2 

2 y di! ^ Hpit' - t")dt" ^Jok%t 







(3.86) 
(3.87) 



which we will use below. The integrals over u can be done by closing the integration 
contours in the complex plane and by evaluating the residues. The final answers in 
formulas (|3.85| )- (|3.87| ), which are given after the right arrows give the results in 
the limit t ^ r, t ^ g^^. 



Using equations (|3.81| ) and (|3.84|) , we easily obtain 



rd 



Jpk '^Ofc (^ + ^rd)^ ~ ^rd 

r 2 + r^d) 



+ 



Jok {Ci + - {2n + n 



rd 



2 (2fi + ad)[l + r(21] + ad)]' 



(3.88) 



Substituting this equation into equation ( |3.53D , we see that the ensemble average of 
"""V squared is finite, but it would be infinite if the rotational damping were absent 
(i. e. if flj-d = 0). 
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Chapter 4 



Energy Spectrum of Random 
Magnetic Fields 



In this chapter we use equations (|3.80| )- (p.83|) for the statistics of turbulent velocities 



to find the evolution of the magnetic field energy and of the magnetic energy spectrum 
in the magnetized turbulent dynamo theory. The principal results of this chapter (and 
of the thesis) are given by the following equations. The growth of the total magnetic 
energy is given by equations ( [4.16|) -( ^.20| ). The evolution of the magnetic energy 



spectrum is given by equations (|4.56|) - ([4.63|) . The evolution of the magnetic energy 
spectrum on small (subviscous) scales is given by equations ([4.88|)-(p:.9g|). 



4.1 The Growth of the Total Magnetic Energy 

We start with calculation of the magnetic field energy growth because it is of the 

greatest interest. The volume averaged and ensemble averaged [] magnetic energy per 
^ Ensemble averaged over the ensemble of turbulent forces. 
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unit mass is 



£ 



def 



1 
L3 



L/2 



-L/2 



k=0 



(4.1) 



where p is the plasma density, and -B^k=o is the k = Fourier coefficient of the 
magnetic field strength squared, B^, see equation ( |A.5| ) in Appendix 0. To find S{t), 
it is convenient to introduce the following symmetric tensor 



Bai3 == BryBfi = B^ b, 



'afi , 



(4.2) 



where hap is given by equation ( p.38| ). The differential equation for Baf^ can easily be 



derived from equation (|2.39|) 



dtB, 



t-Dal3 



BadtBp + BgdtBa = Va.-rBgy + V^.-y-B, 



' f3 ,'y a"/ 



(4.3) 



Now, we solve this equation and equation ( p. 401 ) by making use of the quasi-linear 
expansion procedure, described in Section |3.1| . First, we write the expansion for B^ 



up to the second order, and for Bap up to the first order |, 



B\t) 



''B^ + 'B\t) + ''B\t), 



^B 



^Bap it) . 



(4.4) 
(4.5) 



Here, at zero time 5^(0) = "E^, ^5^(0) = 0, ^5^(0) = 0, 5^/^(0) = and ^5„^(0) = 
0. Second, we substitute these expansion formulas into equations (|2.40|) and (|4.3|) 



^ Note, that below the magnetic field strength squared B'^ is expanded as a whole. Thus, for 
example, the first order quantity in equation (4.4) is not equal to ■'■B • ■'■B, where -""B is the first 
order expansion term in equation (3.1), the later is of the second order. Of course, is equal to 
2(°BiB). 
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We find that the first order equations are 

dt^B^ = 2V,/S,^- (4.6) 

dt^Bap = ^Va,'r^Bi3^ + '^ViS^j^Ba-y — i^V ^Ba(i),^ , (4.7) 

and the second order equation for '^B'^{t) is 

d^B^ = 2 'V^/B^^ + 2 'V^/B^p - ( V/52),^ - (V/52),^ . (4.8) 

Here, in the last two equations we transform the last (convective) terms by making 
use of the plasma incompressibihty condition, V^^a = 0. 

Next, we first integrate first order equation (|4.6| ) in time (with the zero initial 
conditions) and ensemble average the result. Using equation ( |3.8CI|) , we obviously 
obtain 



'B\t)) = 0. 



(4.9) 



Second, we integrate equation (|4.7] ) in time (with the zero initial conditions), and 
then Fourier transform the result in space, r ^ k. We have 



a-f 



K E '^k'a(t')°^k"/37 + ^; E 'v^k'/j(t')°^k"c 
k' k' 

k"=k-k' k"=k-k' 

- ^ ^Vliij{t') °i?k"a/3 dt' 

k' 

k"= k-k' 

t 

K{Sar%, + Spr'b^,)-K%J f ^ 'V^>r{t'fB\. dt' . (4.10) 



Here, to obtain the last line of this equation, we use equation (|4.2|) and we assume 
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that the magnetic field fines can be considered as initiafiy straigfit, '^bajs = const, 
(tfiis is our first working hypothesis, see Section p.2| ). Third, we integrate the second 
order equation ( [4 .81) in time (with the zero initial conditions) and ensemble average 
the result. Using equation ( p.83| ), we obtain 



V(t) 



dt'. 



(4.11) 



Fourth, we Fourier transform this result in space, r — k, and set k to zero. The second 
(convective) term in the brackets [...] in equation ([4.1 1|) gives zero contribution to the 
final result 0, and we have 



(4.12) 



Fifth, we substitute equation ( [4.1(J ) into this last equation, and use equation ( 3.81| ). 
We obtain 



K=o(t)) = 2 Vk=o / dt' / Y.ix'eCVUt'fV^Ut"))dt" 







2Vk=o^/i^A;2 dt' / \HF{t' -t"]VL + VL,A,VL 



'rdj 




J J.// 



+ Hpit' - 1"; n + ad, 2n + ad) 



dt" 



(4.13) 



^ This is not surprising because the convective term can only redistribute magnetic energy in 
space, but can not change the volume averaged energy. 
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where 



7 



(4.14) 



Here, we also use equation (|3.86|) in the hmit t^r. 

Now, we can finally obtain the differential equation for the averaged magnetic 



energy 8. Following Kulsrud and Anderson ||2^, we choose t small enough for the 
quasi-linear expansion to be valid, but large enough for the limit t ^ r to be satisfied. 
As a result, using equations (|4.4|), ( [4. 91) and (|4.13| ), we obtain 



k=0 



(4.15) 



and using equation ( [4.1| ), we finally obtain 



dS 



7 
2Q 



TT 



27r 



J k^Jokdk //i' (l 



^ = 90^/i2(l-/x2) 



1 + ( 1+ 



djJ, . 



(4.16) 
(4.17) 

(4.18) 
(4.19) 



Here, we replace the summation over k in equation ( 4.14|) by integration over k, 
and then use ci^k = 27ik'^ dkdfi. We also use equations ( p.371 ), (|3.38|) and (|3.64| ). 
Equation ( |4.16|) coincides with formula ( 2.16|) obtained by Kulsrud and Anderson in 
the case of the kinematic turbulent dynamo. However, the growth rate 7, given by 
equation ( [4. 17] ) in the case of the magnetized turbulent dynamo, is different from 
the growth rate 70, given by equation ( ^.17|) in the case of the kinematic dynamo. 
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Let us now compare these two magnetic energy growth rates. First, if we consider 
the kinematic turbulent dynamo, then in our formulas we need to take the limit 
— > 0, i7 ^ 0, or alternatively, the limit fird — ^ cxo 0. In this case, after averaging 
over fi, the growth rate ( |4.17D [see also (|4.14| )] reduces to the growth rate (|2.17| ), 



as one can expect. Second, let us use formula (|3.64|) to make an estimates of the 



magnetized dynamo growth rate 7 and of the ratio 7/70- The calculations are given 
in Appendix |C|, the result is 

7 ^ 80 i^-j^j -f, (4.20) 

7 We - (^) (l^) (4-22) 
^ 10, (4.23) 



where to obtain formulas ( |4.21| ) and ( [4.22| ) , we use typical parameters in a protogalaxy 



given in Table |0| (here ^ is the ratio of the total mass M to the baryon mass, L is 
the protogalaxy size). Thus, our prediction is that the magnetic energy growth rate 
in the magnetized dynamo theory is up to ten time larger than that in the kinematic 
dynamo theory. Note, that two different effects contribute to this difference. First, 
the effective viscosity in the magnetized dynamo theory is smaller than the molecular 
viscosity, u^s = (l/5)z/ < z/, this effect makes the magnetic energy growth rate larger 
by a factor of square root of five (this factor was included in Kulsrud et al. ||28|| ) . 
The rest of the contribution to the difference between the growth rates comes from 



4 



In this limit the "driving" terms 3i^(°6a/3/ii/t^/i,i'),/3 and (1/5)i^AJ7q in equation ( 3.67 ), which 
in the Fourier k-space result in the two corresponding driving terms proportional to SI and f2 
respectively [e. g. see equation ( [B.3| )], go away. Taking the limit fird ^ 00 is equivalent to an 
infinitely large damping of the back-reaction velocities. As a result, in both these limits the back- 
reaction velocities are zero, and the turbulence is Kolmogorov. 
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the local anisotropy of the turbulent velocities in the strongly magnetized plasma. 

Please note, that the number of e-foldings of the magnetic energy during the 
collapse time, 7 tcoUapse, given by equation ([4.22|) , does not depend on L, and therefore, 
does not depend on the redshift. 



4.2 The Mode Coupling Equation for 



the Magnetic Energy Spectrum 



In this section we derive the mode coupling equation for the evolution of the magnetic 



energy spectrum M(t, k), given by equation ( |2.18| ). Using the quasi-linear expansion 
formula ( |3.1| ) for the magnetic field, we obtain the ensemble averaged square of the 
magnetic field Fourier coefficient, up to the second order. 



iBk(t)r) = |°Bk(t)P+[(^5,.(t))°sL(t) + c.c. 



(4.24) 



We calculate all terms in equation ( [4.24| ) separately, our calculations are similar to 
those of Kulsrud and Anderson |]2^ . 



4.2.1 The (iBka(*)>°B*^ + c.c. term 

We start with the calculation of the (^5kQ,(t))°B]^Q + c.c. term. We integrate the first 
order equation in time (with the zero initial conditions) and ensemble average 

the result. Using equation ( p.80| ), we obviously obtain {^Ba(t)) = 0. Thus, 



(4.25) 
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and the (^Ska)°6ka + c-c. term on the right-hand-side of equation ( ^4.24| ) is zero. 



4.2.2 The {\^B]^{t)\^) term 

Next, we calculate the (|"'"Bk(t)P) term. First, we integrate equation (|3.11|) in time 



(with the zero initial conditions), and then Fourier transform the result in space, 
r ^ k. We have 

t 

^B^^it) = tk^{6^J^s - S^J^s) J E 'V^"a{t')''B^'5dt'. (4.26) 



" k"=k-k' 



Here we also use the fact that the fluid velocity and the magnetic field are diver- 
gence free, and therefore, fc^°Sk'a = and fc"^V^k"a = 0. The complex conjugate of 
formula ( [4.26| ) is 

t 

'Bl^it) = -iK{5^p5r, - 5rp5^,) I E 'vl,.p{t'rBl",dt". (4.27) 

. i'"' 

k'v=k-k"' 

Using these two equations, we obtain 



(|"'"Bk(t)P) = k^Ki^aP^^S^T-q — ^ar^^pT^-yS — ^a-^^pS^TT) + ^a-y^Pr^S-q) 
t t 

1 1 E °^k'5°5;,„^(iVk"„(t')'v^kiv^(t"))rft'rft" 

i_/ 1.111 



^ ' k"=k-k' klv=k-k"' 

t t 

X E JJ%n\'B^.'\^'V^''a{trK''f5in)dt'dt''. (4.28) 



Here, we use (^^k"a^^k'^/3) ^k",ki^! see equation ( p.81[ ), and therefore, k'"^ = k" and 
k'" = k'. We also assume that ^bap = const (our first working hypothesis), and there- 
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fore, "Bk'^'^-fik'r; = "^(SrjT^k'P- Now, according to definition ( p.l8| ) for the magnetic 



energy spectrum M{t, k) and equation ( [4.24| ), in order to obtain the ensemble aver 



aged energy spectrum, {M{t,k)), we need to integrate (|-'-Bk(t)P) over all directions 
of unit vector k. Equation ( |4.28D has four terms on the right-hand-side. Therefore, 
we have four terms for (|-'-Bk(t)P) integrated over k, 

yp(|^Bk(t)Hrf2k = T -T -T" + T'", (4.29) 



where 



L - 



T= k^[-^) J dk'J k'^ fB^,\^d^k' J fi^d^k J I ('V^k".(t')'V^I"«(i")) dt'dt", (4.30) 



2 oo t t 

T'= fc^(^) Jdk'Jk'^ lO^k.lVk'l^t/'k J l%aCV^''a{t'yK''pit''))kdt'dt'', (4.31) 



DO 



T" = T'* = r, (4.32) 

3 oo t t 

T"= fc'(^) jdk'jk'^ |05kf c/^k'l^'k j I LCK''a{t'yK''p{t''))hdt'dt'', (4.33) 



k" = k - k'. (4.34) 

Here, /i = (k • ^h), see equation ( |3.32| ), and we replace the summation over k' by 
the integration over k' (see Appendix which, in turn, is replaced by the double 
integration over k' = |k'| and over k' = k.'/k'. 

[In equations ( |4.29| )-( p.33| ) and below, in order to shorten our notations, we use 



the same notation for a Fourier coefficient, no matter whether it is a discrete or a 
continuous function of k. For example, the function (|-'-Bk(t)P) in equation ( [4.28|) 
is a discrete function of k, while the function (|-'-Bk(t)P) in equation (|4.29|) , strictly 
speaking, is the appropriately defined continuous function of k, (|"'"B(t, k)p), see Ap- 
pendix 0. However, it is convenient to use the same notation for the both functions.] 
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plane 
perpendicular 

to k' 



Figure 4.1: This plot shows relative position of vectors ^h, k, k' and k" = k — k' in 
space for the mode-coupling kernel ( |4.57| ). The k" modes of the turbulence interact 
with the k' modes of the magnetic field to change the energy in the k modes of 
the magnetic field. (In our case of an initially straight magnetic field vector k' is 
perpendicular to °la because the field is divergence free.) 
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Now, refer to Figure ^.11 . Note that vector k' is perpendicular to the zero order 
magnetic field unit vector because the field is divergence free, k'J^Bi^/a = 0. We 
also have k" = k — k', see equation (|4.34| ). Therefore, the following useful equations 



are vahd, see also Figure 



kk' 
k" 
kk" 

1^" 



1 — fi 



112 



COS 6 , 

|k"| = {k^ + k'^ -2kk'coseY/\ 

k-k" _ k - (k-k') _ k-k' cose 
k" ~ 'W' ~ ' 

k • °b = sin 6 cos ({> , 

Ob ■ k" _ °b ■ k _ k 

k" ~ k" ~ F 
k"^ - k^ sin^ 6 cos^ (f k^ + k'^ - 2kk' cos 9 ~ k^ sin^ 9 cos^ (p 

kf^ ~ k"^ 

k''^ - 2kk' cos 9 + P - k^ sin^ 9 + k^ sin^ 9 sin^ (p 

{k' - k cos 9y + k"^ sin^ 9 sin^ Lp 



% • k" 



sin 9 cos if . 



(4.35) 
(4.36) 
(4.37) 
(4.38) 
(4.39) 



(fk = sin 9 d9 dip . 



(4.40) 
(4.41) 



Using these formulas, equation ( |3.81|) , and equation ( p.87|) in the limit t ^ r, it is 
straightforward to calculate the double time integral terms in equations ( [4.30| )-( p:.33| ), 



t t 







1 + 



2Q"' 



t t 







XCV^"a{t'yV;,.^{t"))kfsdt'dt" = Jo,„t|^l + ^j [^^W' 



k'{k' — k cos I 



-2 



sin 9 cos ip 



t t 



n" 



kaCV^,,^{t'yV^„^{t"))kfsdt'dt" = J,,„t[l + — 







rd. 



(4.42) 



(4.43) 
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k sin 9 sin (p 



{k' — k cos Oy + k'^ sin^ 9 sin^ ip 

2n"y^ k'^ 

k^ 



1 



X 



{k' — k cos OY sin^ 6* cos^ (p 
{k' — k cos OY + sin^ ^ sin^ ip 



(4.44) 



Here, Vt" and i^J'j depend on k" and on /i"^, see equations (|3.37| ), ( p. 3^ ) and (|3.64| ). 
In turn, k" and on /i" are functions of fc. A;', 6* and y?, given by equations ( [4.36|) 



and (|09| ). 

Now, we substitute formulas (|]3|), ([4.41|) and (|4.42|) -( ^^ into equations (|4.3q ), 
( [4.31| ) and ( |4.33| ). The factors that we obtain in these equations after integration 
over d?]^ = sin 6 dO dip depend only on k and A;', so they can be exchanged with the 
integrations over c/^k'. Combining the results together in equation ( |4.29| ), we obtain 



k'{\^B^{t)\')d'\i 



dk'}Ci{k,k') / jo^kf rf'k' 



(4.45) 



where 



)Ci{k,k') = tk^ 



J de sm'eJok" J dip (i + ^J 

n n V I'd / 



X 



L 



k'^ + 2k{k - k' cos 9) cos^ ip 
-^2 



-2 



(fc'-fccose)2 + (F-fc'2)sin'^sin> ^ 
[k' ~ k cos dY + k'^ sm^ 9 s\.n^ ip ^ 



(4.46) 



and k" is given by equation ([4.36|) . 
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4.2.3 The (2Bka(*)>°B^^ + c.c. term 

Next, we calculate the {'^Bi^ait))^Bi^^+c.c. term of expansion equation (|4.24]). First, we 



integrate equation ( ^.ISp in time (with the zero initial conditions), and then Fourier 



transform the result in space, r — ^ k. We have 

t 

2Sk,(t)=^fc.(Mrx-M.x) / E ['V^k"/3(t')'^k'x(i') + '^k"/3(t')°^k'J^^t', (4.47) 





^ k"=k-k' 



where we use the divergence free conditions kJ^B^^a = 0, ka^B]^a = 0, ka^V^a = 
and ka'^Vka = 0. Second, we ensemble average his equation. The second term in the 
brackets [...] averages out because of equation (|3.83| ). Then, we multiply the resulting 



equation by ^B^^^, add the complex conjugate, and use formula ( [4.26|) for ^-Bk'x- We 
have 

t 

CB^„)'BI^ + C.C. = lk^{5^p5r^~5rM J Yl CV^k"/3(t')'^k'x(0)°^I^^'^' + C.C. 



n 

^ k"=k-k' 



X 



t t' 

Y: Y K fdt'f dt" CV^..^{t") 'V^.p{t')) 'Bl,'B^"'5 + c.c. 



W k' 
k"=k-k' kiv=k'-k 



t t' 

X Y KJ dt'J dt" (^K_k"a(t") 'K"p{t')) + C.C. (4.48) 

k-'k-k' 

Here, we use (^V^k'^a^^k"/?) ^k^^.-k"; see equation ( p.81| ), and therefore, k'"^ = — k" 
and k'" = k. We also use the field divergence free condition, fca^Bka = 0, this is why 
we have only three terms at the end. Third, we change the summation over k' to 
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summation over k" in equation ( [4.48 



). We have 



t t' 

xY^ik,- k'^) J dt'J dt" (^V^_k"«(t") 'V^np(t')) + c.c. 





t t' 



- 2k^kp |°BkP E / dt'j dt" C^-k"a(t") 'K"p{t')) 





Oh |05_ |2 



k" 



X 

k" 





t t' 



- 2e |°BkP T.fdt'f dt" LCV.^.o.{t") 'V^.p{t'))kp 

k" 

t t' 

+ 2k i°Ekp E °^7^;' / dt'J dt" KCv.^>un 'v^"p{t'))% 





t t' 



k" 

, 3 oo 

_2fcTBk|'(^) / rf'k" 

— oo 

t t' 

X l^t' 1 rft" L('V^-k"a(t") 'V^"p{t'))kf, . (4.49) 







Here, we again used the divergence free conditions on the field and on the fiuid 
velocities, and formula °Bk,,°Bk5 — °^r)5|°5kP (°&r?<5 = const, according to our first 



working hypothesis). Two terms on the seventh and eighth lines of equation ( [4.49|) 
cancel each other because of the symmetry of tensor (^V_k"ci ^V^k"/?) with respect to 
exchange a ^ j3 and k" ^ — k", see equation ( |3.81| ). We change the summation over 
k" to the integration over k" in the last line of equation ( [4.491 ). 



Now, we integrate (^B\^j^)^B^^^ + c.c, given by equation ( [4.49] ), over all directions 
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of unit vector k. We have 



271 



t t' 

I I t /1t7 U 



X jdt' jdt" K{^V.^.^{t"yV^.p{t'))k^. (4.50) 





Let us refer to Figure ^.21 . Note that vector k is perpendicular to the zero order 
magnetic field unit vector °b because the field is divergence free, /ca^Ska = 0. We 
also have 



K ■ k = cos 6*, 



/i = k ■ b = sin ^ cos 
rf^k" = k"Uk" sine ded^. 



(4.51) 
(4.52) 
(4.53) 



Using the first two of these equations, equation ( p.81| ), and equation ( p.86| ) in the 
limit t r, we calculate the double time integral in equation ( [4.50| ), 



t t' 



dt' jdt" k^{'v.^.^{t"yv^.p{t'))kp 





1 / VL" \ 



X 1- 



X 



1-1 + 



29." 



-2 



COS^ 9 COS^ 



cos^ 6 cos^ if + sin Lp 



(4.54) 



Here, fi", Vt" and 9'^^ depend on k" and on /i"^, see equations (|43^ ), ( p37|) , ( p^Sj) 
and (|3.64| ). Finally, substituting equation ( [4.54| ) into equation ( [4. 501) , and using 
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equation ( |4.53| ), we obtain 



2tt 







rVofc" dk" J sin^ede J d<f 



1-1 + 



2Q" 



X 



cos^ 9 cos^ ip 



cos^ 9 cos^ + sin 



(4.55) 



4.2.4 Collecting the terms together 



Next, we substitute equations ([4.24|) , ( [4.45| ) and ( [4.55|) into equation (|2.18|) for the 
magnetic energy spectrum M(t, k). We also make use of equation ( |4.25|) . We choose 
t small enough for the quasi-linear expansion to be valid, so that dt{M{t, k)) = 
[M{t,k) — M{0,k)]/t. As a result, we finally obtain the mode coupling equation for 
the magnetic energy spectrum. 



dM 
~dt 



j K{k, k')M{t, k') dk' - 2^ k'^M{t, k). 



(4.56) 



Here the mode coupling kernel K{k, k') is 



K{k,k') 



k'^ + 2k{k- k' cos 9) cos^ (p 



27r 



X 



1-1 



2Q" 
O" 



k" 



{k' -kcos9y + {e-k'^)sm'9sm'ip ^ 
{k' — k cos 9)'^ + k"^ sin^ 9 sin^ (p ^ ' 

(/t=' + A;'2-2A;A;'cos^)l/^ 



(4.57) 
(4.58) 
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k'^ {k' - k cos &f + k^ sin^ 6^ sin^ 



sin^ 6' cos^ 



(4.59) 
(4.60) 



and the turbulent diffusion constant 



47r 



1 r 3 °° 277/ - , 

n n n V 



rd 






2fi 



COS^ COS^ 



cos^ Q cos^ + sin ip 



sin 6' cos V9 (1 — sin 6* cos </)). 



(4.61) 

(4.62) 
(4.63) 



Tlie function Jq^ is given by equation (|2.11|) . To derive tlie above formulas, we use 
equations ( |3.37|) , (p.38|) and (p.64|) . To obtain equation ( [4.601 ) we use equations ( [4.39| ) 
and ([4.40|) . To obtain equation ( ^4.63|) , we use equation ( [4.52| ). In the last three equa- 
tions for the turbulent diffusion constant we replace k" by k'" in order to distinguish 
it from k" in the equations for the coupling kernel. 

If we consider the kinematic turbulent dynamo, then we need to take the limit 
n ^ 0, i7 — > 0, or alternatively, the limit i^rd ~^ oo (see the footnote ^ on page |5^ ). 
In this case, as one might expect, after integrating over y), equation ( f4.57| ) reduces to 
equation ( p.22|) , and after integrating over both Q and equation (|4.61|) reduces to 
equation (|2.24| ). 
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4.3 The Magnetic Energy Spectrum 
on Subviscous Scales 



Equations ( |4.56| )-( |4.63| ), which give the evolution of the magnetic energy spectrum 
in the magnetized turbulent dynamo theory, are the principal result of this thesis. 
However, these equations are rather complicated for easy interpretation. In this 
section we limit ourselves to the evolution of the magnetic energy spectrum on small 
subviscous scales. In this limit, k ^ k^, and the integro-differential equation for the 
magnetic spectrum evolution simplifies to an ordinary differential equation. At the 
same time, the evolution of the magnetic energy spectrum on subviscous scales is of 
great interest for the origin of cosmic magnetic fields (see Chapter |^). 

Let us refer to equation ( [4.57| ) for the mode coupling kernel K{k, k'). The function 
Jofc" cuts off at the viscous wave number k^, [see equation ( |2.11| )]. Therefore, in the 
large-fc limit, k 3> k^, we have k" ~ |A; — k'\ <^ k,k', and can expand the kernel 
K{k,k'). However, the simplest way of calculations is to introduce an arbitrary 
function of k, F{k), which varies slowly in the region k ^ k^, and vanishes outside of 
this region To derive the mode coupling equation on small (subviscous) scales, 
we calculate the following integral 



1 rL 



Airp 



dt 



t=o 



Anp \27r/ 



7T- m 



(|B,(t)p)-|B,(0)p 
t 



F{k) il'B^it)]') d'k 



4:7ip\27lJ t 

oo 

+ / F{k)\CBUt))'Bla + c.c. 



d^k 



(4.64) 
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To obtain the first hne of this equation, we use equation ( p.l8| ). To obtain the 
second line, we replace the time derivative by the time finite difference, assuming 
that t is small, and our quasi-linear expansion is valid. To obtain the final result in 
equation ( [4.64|) [the third and the fourth lines], we use equations ( |4.24|) and (|4.25| ). 
Next, we use equations ( [4.29|) -( p:.34|) and equation ( [4.50| ) to obtain the term in the 
brackets {...} in equation ( [4.64| ), 



F{k){\'B^{t)\')d'k+ / F{k)\CBUt))'Ka + c-C- 



where 



F{k)dk km^B^{t)\^)d^k+ F{k)dk e [fSk„(t))°5L + c.c. 



27r/ 



(4.65) 



T' 



T" 



oo oo t t 

,2l2z?/L\ |0d_ .|2 ^3i,/j3u /"/"/I 



fi''k''F{k) |"5kf rf^k'd^k J J {'Vk"a{t'yV^„^{t")) dt'dt" 



oo t t 

'5kf rf'k' / fi"^k"^F{k) d^k" f [ {^V^u^{t'Yvl„^{t")) dt'dt", (4.66) 



— oo 
oo oo 



-oo 







t t 



l2kF{k) |°5kf c/^k'rf^k J J %CVk"a{t'yK"(iit"))ki3 dt'dt" 

— oo — oo 

oo oo t t 

°5kf ti'k' J fi"k"F{k) d^k" J j %^{^V^n^{t'f Vl,,p{t"))k'p dt'dt" , (4.67) 

— oo — oo 

TL (4.68) 



'O 1 

oo oo 



t t 



F(A;) l^^k'prf^k'rf^k jj k^{^V^.^{t'fvl„ pit" ))kp dt'dt" 

-oo — oo 

oo oo t t 

'5kf rf'k' jF{k) d^k" J I k' J^V^.ait'Yvl" pit")) k'p dt'dt", (4.69) 

-oo 



k = k' + k 



oo 



(4.70) 
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and 



oo t t' 

2 J F{k)fB^\^d'k J d'k" Jdt'Jdt" KCV^^.o.{t")'V^"p{t'))kp 

— oo — oo 

oo oo t t 

J F{k) |°fikpt^^k J d^k" J I k^CV.k"a{t") ^V],np{t'))kpdt'dt". (4.71) 



-oo 



Here, to obtain the final results in equations ( [4.66| )-( p:.69| ), we change the integration 
over k in these equations to integration over k" by making use of formula ( |4.7CI| ). 
We use the divergence free condition on the fluid velocities, k'^Vi^na, and therefore, 
ka^Vk"a = k'^Vk"a- Wc also use the obvious formula 



^k = %-k = %-k" = A;"(°b • k") = fx"k" = k" sin 6 cos 



(4.72) 



see Figure To obtain the final result in equation ( [4.71| ), we make use of equa- 
tion ( p.81| ) and of the first equality in equation ( p.87| ). 

We calculate 7^, T^" and T" separately, up to the second order in k" <^ k, k'. 

First, following Kulsrud and Anderson in equations ( |4.66| )-( |4.69| ) we expand 
the slowly varying function F{k) in k" <^ k at point k' up to the second order. We 



have, see Figure 4.3 



k" + k'" + 2(k' ■ k 
k' 2 



1/2 



k^ 



= k' 
'k' ■ k"? 



Ill 



2(k' ■ k") k 
1 + -^^ + - 



1/2 



k'^ 



k 



111 



k' + k"cose+'-^sm^e, (4.73) 
2k' 



F(k) 



dk' 

dF I dF Id'^F 

"-^ - " " ± or . „2 ^.^2 , " ul'2 _„2 Q 



F(k') + —k" cose + ——k"^ sin^ 6 + -^k"^ cos' 6 



2k' dk 



2d¥ 



Second, we calculate given by equation ( |4.66|) . Because is of the second 
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plane 
perpendicular 

to k' 



Figure 4.3: This plot shows relative position of vectors °b, k', k" and k = k' + k" 
in space for equations ( [4.66|) -( p:.70D . 
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order in k'\ we need to keep only the zero order term in expansion ( ^.74| ) for F{k). 
Thus, we have 



t t 



% = J F(A;')l°^kf J ii''^k''^d'\^'jjCV^,,^{t'fvln^{t''))dt'dt 

— oo — oo 

OO oo TT 

= t f F{k') dk' f jo^kf d^k' f k"^jQk" dk" f de sin^ 6 















27r 



X Jd(p cos^ '^[^^Q/r 



1+1 + 



2Q" 



(4.75) 



where we use 



d^k" = k"'' dk" sm e de dip 



(4.76) 



and equations ( ^.42| ) and ( [4.72| ), see Figure |4.3| . Here and below in this section, 
functions H", Q!' and QJ.'^ depend on k" and on /x"^ = sin^ 6 cos^ see equations p.37| ), 
O) and { ^W^ . 

Third, we calculate 7^', given by equation ( ^.67] ). Because ^"k" is of the first order 
in fc", we need to keep only the zero and the first order terms in expansion ( |4.74 ) for 
F{k). We have 



t t 



X = F{kXB^'\'d^^' U"k"d'k" %aCV^''a{t'YVl''0{t''))kUt'dt'' 



-oo 
oo 







t t 



+ / ^fB^,\^d'k' ffi"k"\osed^k"l l%aCVk''a{tfK^ (4.77) 



dk' 







Using equations ( p.81[ ), ( |3.87D , ( [4. 72] ) and Figure we obtain 



t t 



\CK''a{t'YK''Bit'')Kdt'dt 
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-tk'Jok" sin 9 cos 9 cos (/3 ( 1 + ^ j ( 1 + -^1 



(4.78) 



Now, we substitute this equation and equation (|4.76|) into equation ( [4.77D . The first 



term in equation (|4.77| ) vanishes after the integration over 9 because the integrand is 
an odd function of cos^. f\ The second term is nonzero. As a resuh, we have 



r; = -tlk'"^ dk' J k'^ |°5kf f/'k' J k"^Jok" dk" J d9 siir" 9 cos^ 9 





(4.79) 



Fourth, we calculate T^", given by equation ( [4.69|) . We need to keep all terms in 
expansion ( [4.74|) for F{k). We have 



T" 



oo t t 

F(A;')r5k'|'c/='k' U^'k!' f f k'^{W^,n^{t'fVl,,.{t''))k'dt'dt 



-oo 
oo 



-oo 

oo t t 



+ / ^r^kft^'k' I k'' cos9 d^V I j k'^{^V^n^{t'fvl,p{t''))k'pdt'dt'' 



-oo 
oo 







-17 7-1 °° * * 



oo 



+ 



2k' dk 



-^51^|"5kfrf^k' /^''^sin^^c^^k'' // A::,(iVk".(tOVk"/.(^'0)^;5rf^^^^^^ 





t t 



2dk'^ 



Xfrf'k' / k"^ cos^ 9 d'k" j j k'^Cyk.a{t'yv^„f,{t"))k'pdt'dt". (4.80) 





The first term in this equation is equal to minus T", which is given by equation (|4.71| ), 
[to see this, we simply change the notation in equation ( |4.71| ), k ^ k']. To calculate 
the last three terms in equation ( |4.80| ), we again use equations ( p.81| ), ( |3.87| ), ( [4.72[ ) 



^ This is because the Braginskii viscosity is invariant under reflection b ^ — b. See also the 
footnote |l^ on page ^ 
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and Figure 0| to obtain 



t t 



K{'v^„^{t'Yvl,„p{t"))k'pdt'dt" 







tfc'Vofc" sin2 6' 1 



1 - 



1-1 + 



20." 



-2n 



COS^ 9 COS^ (f 



cos^ 6 cos^ + sin 



•. (4.81) 



[Compare this equation with equation ( [4.54|) and Figure [4.3| with Figure |4.2|| . Now, 
we substitute this equation and equation (|4.76| ) into the last three terms of equa- 
tion ( |4.80| ). The second term of equation ( [4.80| ) vanishes after the integration over 
9 because the integrand is an odd function of cos 6*. The last two terms of equa- 
tion ( [4.80| ) are nonzero. As a result, we have 



t r,,dF 



2 J dk' 



2tt , 



+ 1 Ik'^dk' I k'^ |°5kf c/'k' / k"^Jok" dk" I de sin^ 6 



rd / 



1-1 + 



2Q" 



-2 



cos^ 6 cos^ ip 



cos^ 9 cos^ V9 + sin 



+ i [k'^fl. dk' I k'^ |°5kf rf^k'y rVofc" dk" I de sin=^ 0cos2 e 





2J' dk'^ 



27r 



1-1 



2Q" 



cos^ 9 cos^ (p 



cos^ 6 cos^ (fi + sin ip 



(4.82) 



Now, we substitute equations (HI), (|475| ), ( ^^91) and j ^M ) into equation (l465| ) 



The first term of equation (|4.82| ) cancels the T" term in equation ( [4.651) . Then, we 
substitute the result that we get in equation ( [4.65|) into formula ( [4. 641) , and using 



equation (|2.18| ) for M(0, k'), we obtain 



Fik) 



dM 



dt 



dk 



t=o 



2n 



dF d'^F' 
XoF{k') + X,k'^ + X.k"'' ^ 



dk'^ 



M{0,k')dk', (4.83) 
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where 



27r 



An 



Ai 



1+1 



j k"^Jok" dk" J de sin^ 6 j dip cos^ ( ^ + 

\ rd, 

y k"\jQk" dk" j de sin^ 6 J dip il + —j I 2 cos^ ^ cos^ ip + - sin^ ^ 
^ '■'^^ I 



(4.84) 



1-1 



2Q" 



-2 



2 cos^ 6' cos^ if + sin 99 



cos^^cosVh (4.85) 



A2 



k"^Jc,k"dk" / sin^^cos^^ 



27r 



1 + 



-2n 



cos^ 9 cos^ 



cos^ 6 cos^ + sin ip 



■. (4.86) 



Next, we integrate the right-hand-side of equation (^4.83|) by parts over some extent 
in k' and choose F{k'), so that it and its derivative dF/dk' vanish at the end points. 
We have 



dt \2txJ 




dM 



(Ao - Ai + 2A2)M + (4A2 - Ai)fc^ + ^2k^-Q^ 



Fdk. (4.87) 



This equation is vahd for an arbitrary function F. As a resuh, the integrands on 
the left- and right-hand-side of this equation should be equal, and we finally ob- 
tain the mode-coupling equation for the magnetic energy spectrum M(t, k) on small 
(subviscous) scales 



dM 



r 

5 



k^^^ - (Ai - 1)A;^^ + AqM 



dk"^ 



dk 



(4. 
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where 



-f-T 

2 \2-n) 



k'^Jok dk I dO sin^ 6 cos^ 6 



2lT 



X j {l + 



^rd , 



to 



















cos^ 9 cos^ ip 



cos^ 6* cos^ + sin if 



Ai 



^ + rid 



Jk^Jok dk JdO sin^e J dip (l + 



X <( 2 cos^ 6 cos^ '/^ + ^ ^^^^ ^ 



1 



sin^^ 



X |2 + - — — " . , cos^gcosV 
2 cos"' y cos"' (fi + sm y 



3 OO TT 2-K , 



X <( 2 sin^ 6* cos^ ^ ~\ ^^"^ ^ + 





1- (1 



^rd , 



2fi\~2 



X cos d - sin^ d 



1 sin^ 6^ cos^ Q 

2 cos^ 6' cos^ <y9 + sin^ 



cos Lp 



^rd 

2fi 



90 — sin^ 6 cos^ (1 — sin^ 6 cos^ 9?) 



(4.89) 



(4.90) 



(4.91) 

(4.92) 
(4.93) 



and function Jofc is given by equation (|2.11| ). Here, we drop double primes, ", and 
use equations (^:84D-(CT) and equations (|337D , j ^M ), (IH- We calculate T and 
dimensionless number Aq and Ai numerically in Appendix (^, and obtain 



r 


^ 100 


Ai 


^ 2, 


Ao 


^ 5. 



(4.94) 
(4.95) 
(4.96) 
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Therefore, 



dM 



20 



L 



k^^^ - k-^ + 5M 



dk^ 



dk 



(4.97) 



It is interesting to compare it with equation ( p.25| ), obtained by Kulsrud and Anderson 



in the kinematic dynamo case [27 



Now, assume that M(t, kj-d) is known as a function of time at some reference wave 
number k = kj-^f, then the solution of ( [4.88[ ) is 



M{t, k) = f M{t\ k,,{)G{k/k,,i, t - t') dt' 

J —oo 



(4.98) 



where the Green's function G{k, t) is 



G{k,t) 



/gAi/2 Ij^fe (r/5)(Ao-A2/4)t -51n2fe/4rt 

5 ^in^ 



An J ri/2t3/2 



e(4r/5)tg 



-51n2fc/4r< 



(4.99) 



We derive these two equations in Appendix 0, and use equations ([4.95|) and ( |4.96| ) for 
Ai and Aq. We see that a "signal" M{t, /cref), at zero time, will increase exponentially 
as e'-^/^-*'"* and will extend down to the scale fcpeak ~ e*^*^/^-""* fcref , where fcpeak is the peak 
of function kG{k, t), (of course, the field scale can not become less than the resistivity 
scale). As a result, in the magnetized dynamo theory the magnetic energy tends to 
quickly propagate to very small subviscous scales, the same way as it does in the 
kinematic dynamo theory, (this propagation is checked by the resistivity). However, 
the tail of the magnetic energy spectrum on /cj-ef ^k < kp^ak scales increases with the 
wavenumber as oc A; instead of oc k^^'^ in the kinematic theory. Note, that according 
to equations (|4.2CI|) and (|4.94|) , the growth rate of the Green's function, (4/5)r, is 
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approximately equal to a half of the growth rate of the total magnetic energy, 27. 

In the end of this section, let us consider the kinematic turbulent dynamo. In this 
case we need to take the limit 0, ^ 0, or alternatively, the limit fird ^ cxd in 
equations (|4.88|) -( p:.9lD [see the footnote § on page |6^. After integrating over and 
9, we have F = 70, where 70 is the Kulsrud-Anderson magnetic energy growth rate, 
given by equation ( p.l7|) , Ai = 3 and Aq = 6. As a result, as one might expect, in the 
kinematic dynamo case equation (^4.88|) reduces to equation ( p.25| ), and the Green's 
function ( |4.99| ) reduces to formula ( |2.27|) . 
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Chapter 5 



Discussion and Conclusions 



As we discussed in the introductory section, the origin of galactic and extragalactic 
magnetic fields is one of the key questions in astrophysics. This thesis is devoted 
to this question. There are two prevailing theories for the origin of cosmic magnetic 
fields. First, the galactic turbulent dynamo theory, which states that the magnetic 
fields have been primarily amplified in differentially rotating galactic disks after the 
galaxies had been formed. Second, the primordial turbulent dynamo theory, which 
states that the fields have primarily been produced in protogalaxies, undergoing grav- 
itational collapse. It seems that both observational and theoretical results favor the 
second, the primordial dynamo theory. 

In calculations of the magnetic field evolution the previous dynamo theories as- 
sumed the regular isotropic viscosity for the turbulent plasma motions, and this is 
justified because of neutrals. However, in protogalaxies the temperature is so high, 
that there are no neutrals, and the viscosity is dominated by ions. Therefore, in 
a protogalaxy, as the magnetic field strength grows in time because of the dynamo 
inductive action, starting from its initial seed value, the plasma quickly becomes 
strongly magnetized, the viscosity becomes the Braginskii tensor viscosity, and the 
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turbulent motions on the viscous scales become strongly altered from the isotropic 
case. As a result, the turbulent dynamo becomes the magnetized turbulent dynamo. 

In this thesis we have developed a theoretical basis for the magnetized turbulent 
dynamo, which operates in protogalaxies. The results of the kinematic dynamo theory 
already seem to support the primordial (protogalactic) dynamo origin of cosmic mag- 
netic fields p8|. The results that we have obtained for the magnetized dynamo, further 
support this primordial origin theory. This is because the number of e-foldings of the 
total magnetic energy during the collapse of a protogalaxy, given by equation ( |4.22| ), 
is as much as ten time larger than that in the kinematic dynamo theory. Therefore, 
the number of e-foldings in the magnetized dynamo is more than large enough for 
the magnetic fields in protogalaxies to grow from their seed value, provided by the 
Biermann battery, up to the field-turbulence energy equipartition value. The number 
of e-foldings of the magnetic energy on the viscous scale, which is equal to the growth 
rate (4/5)r of the Green's function ( [4.99|) , is less by one half, but it is still sufficiently 
large Q. 

Another of our predictions is that the tail of the magnetic energy spectrum on 
the small subviscous scales increases with the wavenumber as oc ^ [see the Green's 
function (|4.99|) ], instead of oc A;^/^ in the kinematic theory. Therefore, in the magne- 
tized dynamo the magnetic energy is slightly shifted to larger scales as compared to 
the kinematic dynamo case. 

Although the important results of this thesis are convincing, we left out several 

issues in our theory. Therefore, let us itemize and briefly discuss the possibilities of 

Of course, our results (|4l7|)-C23|) for the magnetic energ y grow th rate are sensitive to the value 
of the physical parameter fij-^, which is estimated in equation (3.64). We also left out the finite time 
correlation effects (see discussion below). Therefore, our result for the number of magnetic energy e- 
foldings should be viewed as an estimate, valid within a factor of order two. However, it is important 
that the number of e-foldings that we found in the magnetized dynamo theory is large and is clearly 
larger than that in the kinematic dynamo theory. 
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further research on the magnetized turbulent dynamos. 

• Calculations of magnetic field curvature. 

In our calculations in this thesis we assumed that as far as we are interested in 
the evolution of magnetic energy, we can consider the magnetic field lines to be 
initially straight. This our first working hypothesis. It relies on the assumption 
that in the case of the magnetized turbulent dynamo the magnetic field has 
a folding structure similar to the one that exists in the case of the kinematic 
turbulent dynamo 0. Unfortunately, the calculation methods employed in 
this thesis are not adequate to justify this our hypothesis because of complica- 
tions that arise when one calculates the statistics of the field curvature. It would 
be interesting either to expand our theory to the case of curved magnetic field 
lines, or to carry out numerical simulations in order to check the field folding 
structure in the magnetized dynamo theory. 

• Effective rotational damping of turbulent velocities. 

Another hypothesis that we made, but did not proved, in this thesis was the 
inclusion of the effective rotational damping into our equations as a linear damp- 
ing term. Basically, the rotational damping rate fird, which is in front of this 
linear damping term [see equation ( p.65| )], enters our final results as a physical 
parameter. We obtained only an estimate of this parameter, so our final nu- 
merical results depend on its value. As we said above, the rotational damping 
is associated with non-linear coupling of velocity modes and is very important. 
It restricts the unlimited growth of velocity modes which are perpendicular to 
the magnetic field lines and are undamped by the Braginskii viscous forces. It 
is important to further theoretically study the mechanisms which stop the un- 
limited growth of the undamped velocities, and to carry out MHD numerical 
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simulations, including the Braginskii viscosity. 

• Finite time correlation effects. 

In this thesis we used the quasi-linear (up to the second order) expansion proce- 
dure in time to solve the MHD equations. We assumed that we can choose our 
time expansion parameter t small enough for the quasi-linear expansion to be 
valid, but large enough for it to be much larger than the turbulent eddy decor- 
relation time. This is equivalent to the assumption that the turbulent velocities 
are (5-correlated in time. It is known from the kinematic turbulent dynamo the- 
ory that the effects of finite velocity correlation time can decrease the magnetic 
energy growth rate by a factor of order two [0. Thus, including the finite 
time correlation effects into the magnetized dynamo theory is important but 
not vital, since our calculation predicts a very large number of e-foldings of the 
magnetic field strength. 

• Energy equipartition and the inverse cascade. 

Finally, the Green's function solution of the mode-coupling equation for the 
magnetic energy spectrum on subviscous scales indicates that in the magnetized 
dynamo theory the magnetic energy tends to quickly propagate to very small 
subviscous scales, similar to the kinematic dynamo case. On the other hand, 
the observed cosmic fields have rather large correlation lengths. Therefore, the 
magnetic field lines must be unwrapped on small scales by the Lorentz tension 
forces, while the field energy is transferred and amplified on larger scales during 
the inverse cascade stage. This stage happens when the magnetic field energy 
is comparable to the turbulent kinetic energy, i. e. when there is the energy 
equipartition between the field and the turbulence. Thus, calculations and/or 
numerical simulations of the inverse cascade with the Braginskii viscosity are 
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extremely important for the theory of the origin of cosmic magnetic fields. In 
the end of this chapter let us discuss the importance of the Braginskii viscosity 
for the inverse cascade in more details. 

As we said in the first introductory chapter of this thesis, the theory of the inverse 
cascade in a plasma with the regular isotropic viscosity has a difficulty of unwrapping 
of the small-scale magnetic field lines. This difficulty can be understood as follows [|12 |. 



The equation on the turbulent velocities V in the plasma with the isotropic viscosity, 
including the Lorentz forces, is 

dtVa = -P. + /„ + z/A\4 + -^(B-V)5,-(V- V)V;. (5.1) 

47rp 

[Compare this equation to equation ( |2.31| )]. The (l/47rp)(B • V)Ba term is the mag- 
netic tension force, normalized to the plasma density p. The magnetic pressure term 
is included into the hydrodynamic pressure term, —P^a- The uAVa term is the viscous 
term. We can estimate the velocity of the magnetic field lines unwrapping, Kinwrap, by 
Fourier transforming equation (|5.1j ) in space, r ^ k, and then balancing the viscous 
and the inertial forces against the magnetic tension force. For the isotropic case, the 
viscosity dominates on small scales, and we have 



V„ ~ ^%^«^A, (5.3) 



V 



where Va is the Alfven speed. The Alfven speed Va ~ I'k^ at the time of the field- 
turbulence energy equipartition on the viscous scale, and Va < I'k^ before the equipar- 
tition. If the field lines have a folding pattern (and they do in the kinematic dynamo 
theory), then A; ^ /cn ~ k^, see Figure and therefore, the unwrapping veloc- 
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ity ( p.3| ) is small compared to the Alfven speed, even at the equipartition. In other 
words, since Va ^ V at the energy equipartition, then k\\VA ~ 7 (7 is the magnetic 
energy growth rate), and the unwrapping rate /cyVunwrap is much smaller than 7. This 
means that the field continues to grow on the viscous and subviscous scales even at 
the equipartition. 

However, in the case of the magnetized turbulent dynamo, the viscosity term in 
the above equation is modified and the anti-unwrapping argument does not apply. 
Indeed, the field unwrapping velocity is parallel and varies perpendicularly to the 
magnetic field lines. The large velocity gradient perpendicular to the magnetic field 
lines, which leads to a large perpendicular stress in the isotropic viscosity case, is 
irrelevant in the case of the Braginskii viscous forces (because the transfer of the ion 
momentum in the perpendicular direction is inhibited). Therefore, in the magnetized 
dynamo theory /i;||Vunwrap ~ A;||Va ~ 7 at the equipartition, and the magnetic field 
strength saturates on the viscous and subviscous scales. As a result, the Braginskii 
viscosity makes the inverse cascade of the magnetic energy more likely, because the 
larger turbulent eddies can not deliver their energy to the magnetic field on the viscous 



and subviscous scales 26 
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Appendix A 

Fourier and Laplace 
Transformations 

Let us consider an arbitrary continuous function f{t,r), where t is time and r is 
position in space. We assume that / is defined inside a spatial box of length L and, 
in the rest of space, is periodic with the box size. Then f{t, r) can be expressed in 
terms of its Fourier components f{t,k) as 

k 

Here the summation is done over discrete values of vector k, so that the x-, y- and 
z- components of k are discrete and are given by equations 

ky{ny) = -j^riy, (A.3) 

271 

kziriz) = —riz, (A.4) 
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where n^;, Uy and are arbitrary integer numbers. The Fourier coefficients /k(^) are 
given by the standard discrete Fourier transformation formula 

L/2 

h{t) = ^ J f{t,v)e-^'^^dh. (A.5) 

-L/2 

Now let us introduce a function /(t,k'), which is a continuous function of the wave 
number k', 

/(t, k') Ut) if A;,(n,) - ttL-i < fc; < + 7rL-\ (A.6) 

kzi^z) — vrL^^ ^ k'^ < k,j.{nz) + tiL^^. 

This continuous function is constant over small three-dimensional cubic volume ele- 
ments {2tt/L)^ in k-space, and it is equal to the values of the discrete function /k(t) at 
the center points of these volume elements [|. As a result, the Fourier series ( |A.1| ) can 
be rewritten, using the integration of this appropriately defined continuous function 
/(t, k) over k, as 

r 3 °° 

fit, r) = ^ Ut) e^*^^ = (-) / fit, k) e^"^"- d'k. (A.7) 

^ -oo 

In this thesis we frequently use the appropriately defined functions that are continuous 
in k, instead of discrete Fourier components. 

We can further Fourier transform functions /k(^) and f{t,k) in time, by making 

use of the standard continuous Fourier transformation t u. For example, for 

^ Note, that the important property /(t, — k) = /*(i,k) is an exact consequence of /_k = /jj 
everywhere in k-space but on the faces of the cubic volume elements. This is not a problem though, 
because the faces occupy an infinitely small three-dimensional volume. 
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function fk{t) we have 



oo 

Mu) = ^ [ k{t)e'-'dt, (A.8) 



27r 

— oo 



1 

= ^ / /kMe-^'^^da;. (A.9) 

— oo 

We can also form the continuous Laplace transformation in time, i — > s, 

oo 

/k(s) = jk{t)e-^'dt, (A.IO) 



^(^) = «/k(«)-/k(t = 0), (A.11) 

c+ioo 

/kW = ^ / /k(«)e*Ms. (A.12) 

c— ioo 

In the last formula all poles of function /k(-s) should be to the left of the integration 
contour chosen in the complex plane. 
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Appendix B 



The Ensemble Average of the 
Second Order Velocities, (^Vka(*)) 



To derive equation ( |3.52| ) for the ensemble averaged second order velocity in our 
case of an initially straight field ^bap = const, we proceed as follows. Let choose a 
system of coordinates in which the initial field '^b is along the x-direction. In this 
case °6o/3 = Sax^iSx, and we rewrite equation ( p.20|) as 



+ 3l^Sax[\uC'^^,,^ + U^,u)]^^ - [^VaUf3 - UjVf3 - ^Va^Vp]^^ , (B.l) 



Here, we use formula ^hap^u = ^ba/s'^bf^^ + ^ba/3^b^^, see definitions ( ^.3(i| )- p.38| ). Now, 
first, we Fourier transform this equation in space, r ^ k. Second, we multiply the 
transformed equation on the left by the tensor 6;^^ = 6^a — k^ka to eliminate the 
pressure ^P, and we also use the fluid incompressibility condition kJ^Va = 0. Third, 
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we ensemble average the resulting equation. As a result, we obtain 



k' 

k"= k-k' 

k' 

k"=k-k' 

'I3^yk''a) + {Ukia^Vki'/s) + {^Vk'a^Vk'i/3) \ . (B.2) 



k' 

k"= k-k' 



Next, in order to calculate the first order quantities ^&ka/3(^) and ^V]ia{t), which en- 
ter equation ( |B.2| ), we proceed as follows. First, we Fourier transform equation ( p.l3| ) 
in space, and we multiply the transformed equation on the left by tensor 6^^. We 
have 



(B.3) 



Integrating this equation with zero initial condition, ^v\^^q = 0, we find 



L fj, 



uut') + 'yUt')\dt' + nj f/k^^jt') dt' , (B.5) 





(B.4) 



where index 7_|_^ is equal to y or z, and the viscous frequencies Q and Q are given by 
equations (p.37|) and ( p.38| ). Second, we use equation (|3.12|) to obtain 
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Now, we Fourier transform this equation and integrate the resulting equation. We 
have 







dt'. 



(B.7) 



Finally, we substitute formulas ( |B.4| ), ( [B.5| ) and ( [B.7| ) into equation ( [B.2| ), and 
carry out the ensemble averagings. Because Uk'Uk" oc 5k", -k'? we find that k = k'+k" 
is equal to zero, k = 0. Thus, all terms in the right-hand-side of equation ( [B.2| ) vanish. 
Because (^V\^^) is initially (at t = 0) is zero, it stays zero in time, (^V]^^{t)) = 0, and 
using equation (|3.22|) , we immediately obtain formula (|3.52|) . 
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Appendix C 



The Calculation of the Magnetic 
Energy Growth Rate, 7 



It is convenient to split the right-hand-side of equation ( [4.17] ) into two terms: 



7 = y + 7", (C.l) 



7' = TT 



— j J k'j,,dkj fiUl + — \ dfi, (C.2) 

-1 ^ 







Using equations ( p.37| ), ( |3.38| ) and ( |3.64| ), we have 
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First, we calculate the first term, 7', given by equation ( |C.2| ). Integrating over /x, and 
using equation ( |2.11| ) for Jofc, we obtain 



00 







fco 

/ , \ 2/3 1 

= ^f^oA:o(^|^j / x-i/^(l + 6xf rfx 

« 5?^f/„,„fM='^!l^!C5M)"'^£5. (c,5) 

7 " " V^o/ 7 V 27rz/ ; L ^ ' 

Here, we substitute a; = k/k,^ and use the fact that the integral is dominated by the 
upper limit. We also use formulas fco = 2tt/L (L is the system size) and fc^/fco = R^^^, 
where R = Uo/k^VeS = ^Uo/kou is the Reynolds number, and Ues = (l/5)z^ is the 
effective viscosity (see Section 

Second, we calculate 7", given by equation ( p.3|) , in a similar way. We have 



2fco^ ^ °^ V ^3 7 / [l + 90(A;/MV(l-/^')f 



2W y V kl 

2.kf 1k^ ^{k/h)-'''' ( 1 + ^-^Xdk j , rf^z 



/ J, \ 2/3 1 1 2 

^f^ofco 7^ / x~'/' (1 + Qx'Ydx [ ^ 2 dij, 

Uo/ / ^ ^ / [l + 90a;V(l-/^')] 

0.367r[/oA;o ( < (C.6) 



Here, we again use equation for Jofc, substitute x = k/k^, and carry out the in- 



tegration numerically (it clear that 7" <^ 7' because 90 is a large number). Combining 
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equations ( |C.5D , ( |C6| ) and ( |C1| ), we obtain equation ( |4.20| ). 

Third, we calculate 70, given by equation ( |2.17| ), in order to estimate the ratio 
7/70- We have 



7o 





Uo 




2ko 


ko 


nUoko 


(k, 
Uo 



\ 2/3 



(C.7) 



In the case of the kinematic dynamo, the Reynolds number is R = Uo/kou = UoL/2'kv. 
Therefore, using ky/kQ = B?/^^ we obtain 



7o = vr(2vr) {--) - 



(C.8) 



Dividing equation ( [4 .201) by this equation, we obtain formula ( |4.23| ). 
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Appendix D 



The Calculation of Coefficients 



r, Ai and A2 in Equation ( | 4.88 | ) 



Before we numerically calculate F, Aq and Ai, given by equations ( |4.89| )-( p:.93| ), let 
us first check equation ( [4 .881) . To do this, we use equation (p.20|) and equation ( [4.88|) 
to obtain the change of the total magnetic energy £ in time, 



d£ 1 /■ aM „ 1 r 

= - / dk = 

dt 2Jdt 2 5 





dk 



1 r 

2 5 

2^8; 



2 + (Ai - 1) + Ac 



Mdk = -(Ai + Ao + 1)^ 
5 



(D.l) 



where the growth rate 7 is given by equation (|4.17| ). To obtain the second line of 
this equation, we integrate the first line by parts. To obtain the last (third) line 
of this equation, we use equations ( |4.9CI| ) and ( [4.9 Equation ( |D.1D coincides with 



equation ( |4.16| ), as one might expect. 

Now, let us refer to equations ( |4.89| )-( P:.93D . Using equation ( ^.111 ) for Jqa:, and 
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changing the integration variables, k ^ x = k/ki,, we obtain 



n 

^rd 



6x' 



— — = 90x sin 6'cos (1 — sin 6* cos (p) 



u J L ' 



T' 



Ai 



An 



^ i^^y J 1/3 J sij,3 Q ^^g2 Q 


X J ll + 




2 



n 

^rd , 



rd 



COS^ 9 COS^ V9 



cos^ 6* cos^ V9 + sin ip 



5 /Svry/^ j 



T'\2 J 



-1/3 dx J de 



X <( 2 COS 6* cos + 2 ^^"^ ^ ~ 







X 2+ - 

V 2 

5 /57r\i/2 r 



sin^^ 



"2 — 1 "^OS^ ^ '^OS^ V9 

2 COS"* 9 COS"* + sin (p I 



2tt 



r V 2 / 



/ X dO sin^ 9 j dip il 



X <( 2 sin^ ^ cos^ ^ ~ 2 ^^^^ ^ ~^ 



X cos^ 9 - sin^ + 







^rd. 



21]\-2^ 

f^rd / 

1 sin^ 9 cos^ 61 

2 cos^ 9 cos^ + sin^ ip ^ 



COS^ (y9 



(D.2) 
(D.3) 
(D.4) 



(D.5) 



(D.6) 



(D.7) 



Here we use formulas k^ = 27r/L (L is the system size) and k„/ko = R^^^, where 
R = Uo/koUeS = ^Uo/k^v is the Reynolds number, and = i^/^)^ is the effective 
viscosity (see Section |2l^). We calculate the triple integrals in equations (p.5|) - ( |D.7| ) 
numerically. We obtain 



r ^ 104, Ai ^ 2.13, Ao ^ 5.21, 



(D.8) 
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which, immediately give us equations ( |4.94| )- ( ^.96| ). 

It is interesting that if i^rd = ^rd(^) [it depends only on k, not on /i^], then in the 
limit i^rd/ J^cskl —>■ (weak rotational damping) we have Ai = 2 and Aq = 5. The 
results given by equation (|D.8D are close to these results. 
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Appendix E 



The Derivation of the Green's 



Function Solution ( 4.99 i ) 



To solve equation ( |4.88| ), we Laplace transform it in time, t s and M{t,k) 



M{s,k). We assume that M{t,k) is zero at t = and that we know M(t, fcref) at 
some reference wave number k = k^^i at t > 0. To simplify our notations we set 
fcref = 1, SO that M(t, fcref) = M{t^ 1). As a result, we have 

^.M^P|^-(A,-1).§^A„M. ,E.l) 

We consider the following solution of this equation: 



M(S, k) = M{S, 1) A;(Ai/2)±V5-/r-(Ao-A?/4)^ (e_2) 



where 



oo 

M(s, 1) = / e-'''M{t', 1) dt' (E.3) 
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is the Laplace coefficient of function M{t, 1). This solution satisfies M\k=i = M{s, 1). 
Now, the minus sign in the exponent in equation ( |E.2| ) should be chosen to satisfy 
the boundary condition at infinity in k, M{s, oo) = 0, [^]. The Laplace inversion of 



M{s, k), given by equations ( [E.2|) and ( [E.3D , is 

c+ioo 

M{t,k) = — f e''M{s,k)ds 

c—ioo 
oo c+ioo 

= 1-Jdt' J M{t', 1) e^(*-*') k^'/^ ^-v/5^/r-(Ao-A?/4) 

c—ioo 

Changing the integration variables 



t'^T = t-t', s ^ 2 = a/Ss/F - (Ao - Af/4), (E.5) 



we obtain 



M(t,A;) = — A;^i/2 /• M(t - r, 1) e^^/^^^^""^?/^)^ rfr / ;2 e^^"/^)"'-^''^'^ (E.6) 
Svri J J 

—oo —ioo 

Here, the z integration can be taken as the imaginary vertical axis in the complex 
2;-plane. Carrying out the z integration, we obtain 

. c ^l/2^Al/2l L * (r/5)(Ao"A2/4)T-51n2fc/4rr 

M(U) = (A) ^^/M(t-r,l)^ -^^ dr. (E.7) 



— oo 



This function satisfies equation ( [4.88 ), but does not satisfy the initial condition equa- 



tion, M(0, k) = 0. It is clear that function 

/ c N 1/2 lAi/2 1„ L i „(r/5)(Ao-A2/4)T-51n2A:/4rr 
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satisfies both equations. Because M{t — r, 1) is zero if r > t, we can extend tlie upper 
integration limit in equation ( [E.8| ) to infinity. Tlierefore, tliis equation is equivalent 
to equations (|4.98| ) and ([4.99|) , after we replace k by k/kref and t — r by t'. 
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